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Chapter One 

INTRODUCTION TO LASER PLASMA HEATING 

In recent years the need has become apparent for a new 
source of power generation. The conventional sources of 
power are either being depleted {coal, petroleum) or else 
require careful disposal of dangerous waste products (nuclear) . 
If controlled thermonuclear reactions can be generated and 
made economically competitive, then there should be no dan- 
gerous waste products . The fuel is deuterium (D) and tritium 
(T) , isotopes of hydrogen. Deuterium is found in water in 
a concentration of one part in six thousand and thus is 
available in limitless quantity. Tritium is not actually 

found in nature but can be produced by using the neutrons 
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generated in the D-T reaction to bombard Li . ' 

The products of D-T and D-D reactions are clean. Ener- 
getic neutrons are produced but these need not lead to radio- 
active wastes if the proper shielding substance is used. One 
of the products of the reaction is a large amount of energy, 
in the million electron volt range. Then even if only a few 
reactions are produced, a large amount of energy can be released. 

• The difficulty in producing a thermonuclear reaction 
is the coulomb force barrier that surrounds any charged 
particle. Coulomb forces are very strong so that the two 
reactants must hit each other at very high speed and pretty 
close to dead center to cause the fusion reaction. Thus very 
high atomic velocities — and hence temperatures — must be present 
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before a significant reaction takes place. For a D-T plasma, 
the temperature must be in the KeV range before significant 
thermonuclear energy is released. For a D-D plasma, the 
requirement is higher by a factor of three. 

There are two examples of power producing thermonuclear 
reactions that have been observed. One is the hydrogen bomb, 
where the necessary temperatures are achieved by a triggering 
nuclear fission reaction. This is an uncontrolled thermonu- 
clear reaction and is not likely as a good source of commercial 
power. The other example is the sun (or any star) which is 
an astronomical thermonuclear reaction. It is in a sense 
controlled since it persists continuously. The goal is to 
produce a small "sun" in a controlled manner so that the 
released thermonuclear energy can be harvested and used for 
other needs. 

In conventional energy release, such as fire, it is a 
simple matter to maintain a continuous reaction. This is 
highly unlikely with a thermonuclear reactor. The contain- 
ment of a continuous reaction at 10 KeV temperatures and huge 
pressures is a frightening thought. For this reason, most 
studies have been directed toward producing a pulsed type 
reactor where each pulse is short enough so that containment 
is possible. A very simple requirement has been calculated 
stating the conditions necessary to achieve gain in a pulsed 
thermonuclear reactor. This is called Lawson's criterion. 3 
Its simplest expression is that the temperature of the reac- 
tants must exceed 10 KeV, and the product of the number 
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density and the confinement time must exceed 10 sec/cm . 

This criterion is based upon the principle that more energy 
must be reaped from the reaction than is required to prepare 
the next reaction. Thus the problem at hand is to heat a 
D-T plasma to a very high temperature, and to maintain it 
for a certain length of time. 

Significant efforts have been made in several countries 
for more than a decade to achieve a controlled thermonuclear 
reaction. This length of time spent proves the difficulty 
in achieving the very difficult conditions of Lawson's cri- 
terion. A number of devices using various principles and 
geometries have been suggested and tried in search of the 
goal. Most of these devices have two things in common; very 
strong electric currents and magnetic fields which are used to 
heat and contain the plasma. Significant progress has been 

made in this direction — especially in the Soviet Union with 
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the Tokamak device* Yet, even this, the most successful 
device to date falls one to two orders of magnitude short of 
both requirements of Lawson's criterion. The Tokamak has 
not been fully exploited yet, but it is still wise to find 
other approaches that may be better. One attractive approach 

that has come to the front in the last three or four years 
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is the use of lasers to hfeat the reactants. ’ ’ 

Ten years ago, the thought of using a laser to heat any- 
thing would have seemed unrealistic since laser power outputs 
were typically very low, on the order of milliwatts. But 
since then, laser technology has advanced at an incredible 
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rate and powerful lasers have become available. Of most inter- 
est to fusion is the development of the so called giant pulse, 
and mode -locked lasers. Giant pulse lasers produce a burst of 
light of microsecond to nanosecond duration. As much as 500 
to 1000 Joules have been produced in this kind of a pulse. 
Kade-locked lasers produce pulse lengths on the order of pico- 
seconds. As much as 80 Joules have been released in this 
kind of a pulse. Thus, the laser is an excellent way to 
rapidly deliver energy. This is essential or else the heated 
plasma will expand and cool before all the energy is delivered. 

Another advantage of the laser is its ability to get 
energy into a magnetic container. Magnetic pressures have 
long been suggested as a means of containing a thermonuclear 
.plasma. Some approaches have used the MHD system {including 
both currents and magnetic fields) to both heat and contain 
the plasma. This requires a great deal of versatility of the 
magnetic field and the currents. However, now the 'heating 
may be left up to the laser so the magnetic field can concen- 
trate on the containment. 

A thermonuclear plasma may also be confined inertially, 
i„e., within a jacket of very heavy atoms. The large mass of 
the atoms would restrict the expansion of the plasma. This 
approach is essentially a small explosion and the products 
must be contained which may be difficult in that the size of 
the explosion may be significant for a positive yield. 

Several difficulties -have been previously found in the 
laser heating of a plasma. The first difficulty is due to 
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the limited confining ability of a magnetic field. This puts 
a limit on the density that can be contained. But lower den- 
sities lead to very long absorption lengths--some even as high 
as many miles. This difficulty has been reduced with the 
advent of powerful gas lasers which are characterized by long 
wave lengths. The long wave lengths lead to shorter absorp- 
tion lengths . 

Another difficulty is that typical estimates of energies 
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required to satisfy the Lawson criterion are 10 to 10 Joules 
per pulse. This is two or more orders of magnitude greater 
than present lasers. 

A third difficulty arises out of the nature of inverse 
brems str ahlung absorption. Inverse bremsstrahlung is a three 
body collision between an ion, an electron and a photon where - 
by the electron gains the energy of the photon. Thus only 
the electrons are heated directly. The deuterium and tritium 
ions 'must be heated by collisions with the hot electrons. 

This equilibration process may well be painfully slow-espec- 
ially for lower densities, and may not even occur before the 

•dynamics dissipates the plasma. 

The purpose of this work is to examine in detail some of 
the phenomena in laser heating to thermonuclear temperatures - 
Much of the work in this field has been confined to estimates. 
Most previous solutions have relied on complex numerical calcu- 
lations which typically leave a dearth in understanding of the 
physical process that produced the results. This study seeks 
simple analytic solutions wherever possible and avoids the 
computer code approach. 
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Standard procedures of analysis begin with the simple 
and once the simple is understood, the more complex is 
attacked. Accordingly, the first study is of a strictly one- 
dimensional plasma heated by a laser. Only stationary heat- 
ing is considered since maximum temperatures are produced if 
the plasma remains stationary throughout the laser pulse. 
Assuming the two temperature continuum equations, analytic 
solutions have been found for all stationary regimes for a 
plasma with arbitrary density profile. 

One among . several heating regimes may appear . One is 
the stationary-frozen-nonconducting regime (where conduction 
will always refer to thermal conduction unless otherwise 
noted) . The heating is sufficiently rapid that the ions 
remain frozen at their initial temperature and little thermal 
conduction occurs. This regime is characterized by two parts. 
For optically thin plasmas, the entire irradiated part of the 
plasma is heated simultaneously, producing a highly nonuniforra 
temperature distribution. An extrema hot spot appears at the 
point of farthest penetration of the radiation (the critical 
density point). For optically thick plasmas, the plasma is 
opaque so that the radiation bores its way into the plasma. 

This phenomena is called a bleaching wave since the plasma 
behind the wave is made transparent by the heating. Conse- 
quently, the temperature profiles present a wave like character, 
the heated region advancing deeper into the plasma with time. 

The stationary-frozen-conducting regime arises for less 
thxk plasmas. In this regime, thermal conduction is very strong 
so that the temperature is nearly uniform. Sven the plasma 
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beyond the light reflection point (critical density) may be 
heated in this case — by thermal conduction. 

The equilibrium counterparts of these two regimes will 
arise for the lower laser energy inputs and for the longer 
pulse lengths, -In these regimes, the electrons and ions 
remain in equilibrium during the heating. The solutions in 
these regimes are essentially the same as those for their 
"frozen" counterparts. 

Of critical importance is the dynamics that follow the 
stationary heating. ' This process appears primarily in two 
regimes. The frozen-conducting regime is characterized by 
cold ions and isothermal plasma motion. This regime is 
unacceptable in that the ions never are heated to thermonu- 
clear temperatures before the motion dissipates the plasma. 

The equilibrium-nonconducting regime is characterized by 
electron-ion equilibrium and no thermal conduction. This is 
the desired regime for the achievement of a thermonuclear 
reaction. 

The classical threshold for the achievement of gain in 
a thermonuclear reactor is the Lawson criterion. Lawson's 
criterion requires the containment of the plasma for a 
certain length of time at a temperature of 10 KeV. Attaining 
this requires fairly large energy depositions, more than is 
available today, but probably in reach of future lasers. • 

Also required are fairly long plasma lengths unless some 
means of -enhancement of the confinement time is found. 

Following the initial step of understanding a simple 
geometry comes the analysis of a more complex one. A geometry 
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which seems appropriate for a magnetically confined, laser heated 
plasma is a long cylindrical shape, the laser being fired into 
the end of the column. The one dimensional results would apply 
to this shape if conditions across the column were uniform. 

This of course will never happen. In fact the density varia- 
tion across the cylinder will cause gradients in the index of 
refraction, so that the beam will be refracted. 

Refraction is a major problem as demonstrated by failures 

0 

in attempts to propagate a laser beam down a plasma column. 

The reason is that the usual density profile in a magnetically 
confined plasma (maximum density in the center) is unfavorable 
—the beam is refracted out of the plasma column. Calculations 
are made of the refraction of the beam in a plasma column with 
a typical density profile. It is found that the beam is essen- 
tially refracted out of the column after moving only a few 
column diameters. Appropriate focusing of the beam helps by 
roughly a factor of three. 

However, if it is possible to generate a plasma column 
with a favorable density profile (a density minimum in the 
center), then the picture is the reverse: the beam is trapped 

perfectly. Not only that, but the effective absorption length 
is enhanced due to the peculiar refractive effects. Some 
unusual features arise in this case. The ray paths are sinus- 
soidal functions, oscillating back and forth across the column. 
Caustics are formed from the locus of points where the rays 
"reflect" back toward the axis. 

The theoretical results on refraction in a favorable 
density profile are useless if it is impossible to create the 
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favorable profile. This type of profile can be created — and in 
several different ways. One way is to catch the plasma in its 
formation stage. Characteristic of the formation of pinch 
devices are collapsing and reflecting shock waves. These 
shock waves generate — for very short times — density profiles 
that are favorable. 

Another approach is to change an unfavorable density pro- 
file to a favorable one by selective heating . This technique 
is called "density tailoring". A laser (called the pre-laser) 
is used for the density tailoring. It is usually advantageous 
to have a shorter wave length for the pre-laser than for the 
main laser pulse. Analysis of the lateral dynamics of a plasma 
column reveals that the motion will fall into one of two regimes. 
There is the conducting regime which is dominated by thermal 
conduction and the nonconducting regime. Several density 
tailoring techniques are suggested and studied. Most of them 
can be used in either regime, but in the conducting regime the 
favorable profiles produced will disappear rapidly when the 
pre-laser pulse ends . 

Oftentimes in the analysis of a particular physical pro- 
blem, the solution is seen to resemble solutions to other 
physical problems. This is a clue to mathematical similarities 
in the governing equations. Occasionally however, the resem- 
blance is not due simply to chance mathematical resemblance, 
but rather to a principle of physics in which both problems 
are rooted. An apparent chance resemblance is seen between 
the bleaching wave and nonlinear thermal conduction phenomena 
(thermal waves) . Further examination reveals that these two 
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phenomena are actually the same but in opposite limits- Both 
belong to a large class of nonlinear diffusion phenomena that 
demonstrate wave like characteristics. The bleaching wave is 
diffusion in the "coherent" limit, and the thermal wave is 
diffusion in the "equilibrium" limit. Two other physical 
phenomena are identified as diffusion in the coherent limit: 
laser irradiation of a two level substance, and the laser 
induced breakdown wave. Numerous examples are presented of 
diffusion in both limits. 

Understanding wave diffusion suggests the possibility of 
laser heating of a very dense plasma by thermal waves. A 
plasma at solid density may have laser energy added at its 
edge by anomalous absorption, and the overdense portion would 
be heated hy a thermal wave. Such high densities would require 
inertial confinement, but the problems of electron-ion equili- 
bration may be overcome by those high densities. Also, the 
high density may permit heating by absorption of a-particles 
from the thermonuclear reaction. 



Chapter Two 

LRSER HEATING OF A ONE -DIMENSIONAL STATIONARY PLASMA 

An important limitation on the maximum temperature 
that can be created in a laser heated plasma arises from the 
plasma motion. During the laser pulse the plasma may begin 
to expand and much of the thermal energy would be converted 
to directed kinetic energy of the plasma motion, as was first 
noted by Basov 5 and Dawson . 6 However, the recent development 
of powerful sub-nanosecond Q-switched lasers and mode-locked 
picosecond lasers makes it possible to overcome the problem 
of expansion during the heating process. For appropriately 
short laser pulses, the plasma can be heated to thermonuclear 
temperatures before significant macroscopic fluid motion begins 
The irradiation of a plasma by^ these short laser pulses then 

becomes a stationary heating problem. 

A study has been made of the stationary heating of a one- 
dimensional plasma with an arbitrary density profile. The 
purpose of the study was to determine the fine details of the 
heating process as well as the overall features. The spatial 
distribution of temperature is one detail of interest which 
will hold the key to the plasma motion that will follow the 
heating process. Another detail of interest is the fraction 
of laser energy that is lost by reflection from the nonuniform 
plasma. 

Several theoretical investigations of laser plasma heat- 
ing have been made; but most models have neglected the fine 
details and concentrated on over-all results, e.g., the tem- 
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uerature achieved. Previous investigations have also used 
time scales in which plasma motion is significant. 

A study by Fader 9 considered the expansion of a spher- 
ically symmetric plasma subject to Q-switched laser radiation. 
The expansion of the plasma was found to be a major limitation 
on the temperatures achieved according to both numerical and 
analytical calculation. This study assumed spatially uniform 
absorption throughout the plasma. 

A study by Kidder 3-0 also considered the motion of a spher- 
ically symmetric plasma but with irradiation by a radially 
convergent light pulse. His numerical study also showed the 
limitation of expansion on the temperatures achieved. 

Dawson, et. al»-, 3-1 studied the over-all absorption char- 
acteristics of a one-dimensional nonuniform overdense plasma. 
Without studying the details of absorption, they found the 
length of plasma necessary to adequately absorb the radiation. 
The critical density is the plasma density at which electro- 
magnetic waves cannot propagate, i.e., the plasma frequency. 
Underdense and overdense refer to densities less than and 
greater than the critical density. 

This study also considers a one-dimensional nonuniform 
plasma and examines the heating process for time scales suffi- 
ciently short that plasma motion can be neglected. The 
results indicate that for a thick underdense region, a 
bleaching wave propagates into the plasma, heating it in 
the process. For a thin underdense region, the heating is 
simultaneous— -with an extreme hot spot appearing at the critical 



13 

density. For a very thin underdense region, thermal conduction 
dominates the heating so that nearly unifirm temperatures appear 
in the underdense region and thermal conduction carries energy 
to the overdense region of the plasma by means of a thermal 
wave. The thermal wave- phenomena has been mentioned by several 
writers (see chapter XV) . The bleaching wave was seen to appear 
in the numerical studies of Kidder , 10 and the analytical studies 

of Rehm „ 1 ^ ' 

Another result is that maximum temperatures or other 
optimizing conditions can be achieved by proper tailoring of 
the initial density profile and plasma size. The idea of tailor 
ing the plasma to optimize the results was first suggested by 
Daiber, et. al ., 7 and later in another scheme by Lubin. . 
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ft.' EQUATIONS AND PARAMETERS 
1. Model 

The plasma is taken to be a one-dimensional nonuniform 
fully ionized gas being irradiated by a laser of wave length 
\ at normal incidence. Figure 1 is a diagram of the config- 
uration. The initial density and temperature profiles are 
arbitrary to the extent that a fully ionized gas is a valid 
assumption. The plasma is assumed to be a dense ideal gas 
mixture of electrons and ions such that there is charge 
neutrality at every point for all times. 

The one-dimensional two- temperature continuum equations 
will be used with the addition of terms accounting for radi- 
ative energy addition. The electromagnetic force will be 
neglected since its effect is small compared with the effect 
of thermal forces for regimes of interest. The effect of 
electromagnetic forces is studied in the first appendix to 
this chapter. Also, a form of the radiative transfer equation 
is added which neglects line, free-bound, and bremsstrahlung 
radiation. The radiative transfer equation is simplified to 
account only for the laser radiation entering the plasma 
(+ x direction) and does not take into account the reflection 
of radiation that would occur if the plasma is overdense. If 
the laser radiation does penetrate to the critical density then 
appropriate calculations must be made to account for the 
reflected light. The absorption coefficient used is that of 
inverse bremsstrahlung . 
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Figure i« Hasma configuration 
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a„ Nondimensionalization Scheme 
The length scale used is R . If the plasma is overdense, 
then x = R will define the point at which the critical den- 
sits’ P c occurs. If the plasma is not overdense at any 
point, then x = R gives the rear edge of the plasma. The 
time scale is the laser pulse length t p since only the heat- 
ing process is considered in this study. Hence, the indepen- 
dent variables, dimensionless distance and time, will always 
be of order one or less. The temperature scale used for both 
electron and ion temperature is T q and is defined as the 
temperature change that would arise if all the energy in the 
laser pulse were added uniformly to a critical density plasma 
of thickness .R . Then, 


f* n ec R 


(2-1) 


where J 0 A is the energy contained in the laser pulse per 
unit area of plasma, k is the Boltzmann constant, and n ec 
is the critical electron number density given by 


ec 


4n 2 c o m e c 0 2 


where e 0 is the permittivity in a vacuum, m e is the mass 
of an electron, c 0 is the speed of light in a vacuum, and 
e- is the charge on an electron. The scaling factor for the 
density is the critical density which is related to the crrti 
cal electron number density by P c = m i n ec ^ * ra i 13 the 
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mass of the ion, and Z is the charge on the ion. The pressure 
scale p Q is the pressure that would arise in a perfect gas at 
temperature T q and density p c . The velocity scale is the 
acoustic speed for a perfect gas at p and p , i.e., . 
a 2 = 5p o /3 Dc . The electron and ion entropies are scaled using 
the constant volume specific heat. The radiation intensity 
scale is the average laser intensity during the pulse, I Q — 
j^/Atp . The absorption coefficient, electron-ion equilibration 
time , and the electron thermal conduction coefficient are 
taken from Spitzer. 14 These are based upon a scale using a 
density p c and a temperature T q 

K 0 = 1.13 x 10" 19 n ec /T^ /2 - • 6.81 x 10 13 Zp c AfT 2/2 m" 1 , 

t e = 3.17 x loV /2 W/n ec Z(l+Z) = 5.26 x 10~ 26 T 3/2 W 2 /p c Z 2 (1+Z) 

690 3 5/2 sec > 

K = 7.40 x 10 T 7 /?, Joule/sec-m-eV . 

e 0 o 

for n ' cm ' 3 . Otherwise, MKS-eV units are used everywhere, 
ec 7 

The Columbic logarithm, InA = 10 everywhere. Henceforth; , the 
ionic charge Z , and the atomic weight W are taken as one and 
two respectively as for deuterium. The radiation intensity scale 
is the average laser intensity during the pulse, I Q = J o /At^ . 

The dimensionless variables are: distance, y = x/R ; 

time, T = t/t p ; density, £ = p/p c ; pressure, n = p/p Q ; 
electron temperature, P e = T e / T 0 < ion temperature, = 

Ti/T ; absorption coefficient, K = k/k 0 ; fluid velocity, 

U = u/a Q ; electron entropy, a e = s e / c v ; ion Enfcr 5>PY> “ 

B-/C ; and intensity, i = l/l Q * 

b. Parameters 

Using the nondimensionalization scheme described, five 
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parameters arise which control the plasma behavior. The first is 


' - e . " 3 o fc p /R (2-2) 

where . e is the ratio of the laser pulse time to the time for 
an electron acoustic wave to traverse the underdense region. 

It is called the dynamics paramter. If e « 1 , then the 
pulse length is so short that essentially no plasma motion 
occurs.- If e & 0(1) , then significant motion occurs before 
the pulse ends. 


<x “ * (2—3) 

a. is the ratio of the thickness of the underdense region to 
the absorption length scale. It is called the thickness para- 
meter since it is a measure of the optical thickness of the 
'irradiated plasma. The absorption length is the reciprocal of 
the absorption coefficient. If a « 1 , then the plasma is 
nearly transparent to the radiation. If the plasma is over- 
dense, then the radiation will be reflected and the reflected 
beam will be negligibly reduce in intensity. The result will be 
to effectively double the absorption coefficient. If a » 1 , 
then the plasma is nearly opaque to the radiation. 


X 


t /t 
P eq 0 


\ is the ratio of the pulse time to the electron-ion equili- 
bration time and is called the equilibration parameter. If 
X « 1 , then the electrons will tend to be heated by inverse 
bremsstrahlung; but the ions will remain nearly frozen during 
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the laser pulse. If X » 1 , then the electron and ion tem- 
peratures will he the same. 


T) 


e.i 


K /pC R 
1 eo c v 


(2-5) 


where r> is the time for significant thermal conduction to 
traverse the underdense region. It is called the thermal con- 
duction parameter. Then, for example, if r) e « 1 , there will 
be very little thermal conduction by the electrons during the 
laser pulse. If q » 1 , the thermal conduction will be 
strong, tending to equalize the temperature of the plasma. 

c. Equations 

The governing equations in dimensionless form are: 


& 
t> T 


+ e 


6£U 

Py 


0 


( 2 - 6 ) 


bcre 

6t 


+ eU 


Sui 

6t 
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by 
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j 
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. i _ 1 

c, 

8 -0. 
e i. 


>2 g7/2 
2/7 5 9 e 

by 

~ ^ re ~ A 
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0 
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+ % 

by 2 - 

+ eU 

boi 

_ x C 

V e i 

' 1 M 

2/7 J 

b 2 0?/ 2 

X 

by 

“ X W2 

6 i 

+ 

ee i 

&y 2 


(2-7) 

(2-8) 

(2-9) 


&i 

by 


+ OKi = 0 


( 2 - 10 ) 


0g /2 (l-C) 1/2 


K = 


( 2 - 11 ) 
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!»=* f V- 

(2-12) 

C 2//3 exp(a e ) 

(2-13) 

C^expfcK) 

(2-14) 


The equation for the absorption coefficient (2-11) breaks do’im 
for densities very near the critical density (C - 1) ♦ There 
is actually an extremely high value of the absorption coeffi- 
cient at that point, but not an infinity. Hence, solutions 
that are found will probably -have a singularity at that point 
and thus cannot be trusted for C very near one. This singu- 
larity will turn out to be integrable so that integrations 
with respect to y over the singularity will describe the 
actual case fairly accurately. 

The initial conditions will be an initial temperature and 
an initial density profile. The boundary conditions will be 
on the intensity; the intensity entering the plasma will be 
d$/dT where $ is the fraction of the laser energy delivered 
up to time T ; 


O 

so that 5=0 at t = 0 and 5=1 at T = 1 (the end of 
the pulse) . 
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2. Preliminary Study of Parameters 
The size of the parameters e , X , r) . will determine 
the character of the plasma heating and the important parts of 
the governing equations. ' One of the factors determining these 
parameters is the size and density of the plasma. For this 
discussion, the plasma is assumed to be overdense {i.e, the 
plasma frequency at the point of highest density is greater 
than the laser frequency) . Furthermore, to simulate a typical 
nonuniform plasma, the density is assumed to rise linearly 
from zero at the edge of the plasma and remain linear at least 
to the critical density point. 

a. Dynamics Parameter 

From (2-2), e = a 0 t^/R . R depends on both the electron 
density gradient and the critical density (which depends on the 
laser frequency) . 


R = 



2 2 

C 4rr e m c 
o e o 

. » 2 


] 




(2-15) 


The acoustic speed scale a depends in the usual way on the 

0 

temperature scale T 0 . A deuterium plasma is assumed in cal- 
culating the acoustic speed. Applying (2-25) T becomes. 



xV e (J 0 /A) 

— am ~2 2 
4rr cm c 
o e o 


2 

e 


] 


Then, e can be written 
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c = 1.04 X 10" 4 (J 0 A) 1/2 X 4 (^n e ) 3/2 t p 

with in Joules per square meter, X Q in microns, 

in electrons per meter to the fourth, and t in seconds. 

« F 

b. Equilibration Parameter 

14 

Using Spitzer's electron-ion equilibration time, with 

(2-4) 

X = 5.76 x 10 66 (J o /AT 3/2 (Vn e f 3/2 x; 8 t p . 


c. Thermal Concuction Parameter 
Using Spitzer’s^ 4 electron thermal conduction coefficient, 
the conduction parameter (2—5) becomes 

= 4.55 x 10~ 148 (J„A) 5/2 (7nJ 9/2 X . 


Then, the pulse time t , laser energy per unit area 

J /&', laser wave length X , and electron density gradient 
0 . o 

7n arise as the key parameters of ’the problem. With a par- 
0 

ticular laser pulse length t and wave length X„ the var- 

P 

ious regimes can be shown graphically as a function of J 0 A 

and . figures 2, 3, 4, and 5 show when e , X , and q e 

—9 —12 

are large or small for the cases t = 10 , 10 sec for 

P 

X o = 1.06(j and t = 10~ 8 , 10 9 sec for X o = 10. 6p. 

The various lines in figures 2-5 separate the different 
regimes from each other. Each regime is characterized by a 
specific combination of large and/or small parameters. Where 
these regimes are distinct, they have been labeled with the 
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Figure 4> Heating regimes for a microsecond 
pulse from a COg gas laser 
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names, “stationary-equilibrium-nonconducting" , "stationary- 
frozen-nonconducting" , "stationary- frozen-conducting" , and 
nonstationary. Conduction will always refer to thermal con- 
duction unless otherwise noted. Since the purpose of this 
chapter is to study stationary heating, the nonstationary 
regime is not divided into specific regimes relating to 
conduction or equilibrium. 

Supposedly in a particular regime, each effect is either 
present or not at all. For example in the stationary-frozen- 
nonconducting regime, one expects that there is no thermal 
conduction, no plasma motion, no ion heating, and only the 
electrons are' heated by the laser. Thus in a particular 
regime, the neglect of certain effects requires that the appro- 
priate parameters be much less than one {or much greater as 
the case may be) rather than just less than one (or greater 
than one) . Hence the approximate solution found will probably 
fail near the lines bounding the various regions in figures 2-5 

It is not always enough for the parameters to be small. 

The frozen ion condition may break down while X is small. 

If the electron temperature is much greater than the ion tem- 
perature, then (which appears in (2-9) ) is not neces- 

sarily small even though X may be small. Also, if the non- 
dimensional temperature at some local' point has a first or 
second derivative not of order one, then the solution may 
not be valid at that local point. For example, in the sta- 
tionary-frozen nonconducting solution, there may be a large 
temperature gradient in some regions of the plasma. Signifi- 
cant thermal conduction may occur at those points even though 
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£be r emainder of the plasma is nonconducting* 
d. thickness Parameter 

At this point it is of interest to calculate the thickness 
parameter - a • The solution may take on a widely differing 
character for different values of a . From (2-3) , 

a = 2.35 x 10 88 ( Jo Ar V2 (Vn e 7 5/ V 10 

In considering the three stationary regimes, assumptions will 
be made — neglecting certain terms depending on the regime. 
However, terms in a will' never be neglected so that a will 
appear in every stationary solution. Thus, in a sense a is 
a free parameter governing the solution, a parameter that may 
be large or small. 

Once having solved for a particular regime, it will be 
necessary to re-examine the various regimes of validity. This 
is bo be expected since the temperature depends on a , and 
thus its order of magnitude may be large or small compared 


to one. 
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B. STATIONARY-FROZEN-NONCONDUCTING plasmas 
1. Equations 

For very rapid energy addition, the plasma may absorb a 
large amount of energy before the random thermal energy begins 
to change to ordered fluid motion. The heating may be rapid 
enough so that the diffusing effect of thermal conductivity 
will be negligible. Also, with very rapid heating, the equili- 
bration processes whereby the heated electrons transfer their 
energy to the ions will not have proceeded to any significant 
extent . 

The problem studied here is the case where all of these 
conditions appear simultaneously. Then, the heating is essen- 
tially a stationary, nonconducting process whereby energy is 
transferred from a beam of light to the electrons in the 
plasma. No energy is transferred between electrons by thermal 
conduction, and no energy is transferred to the ions, either 
by electron-ion equilibration or by ordered plasma motion. 

The requirement for a stationary plasma is that the dyna- 
mics parameter e « 1 . For a nonconducting plasma the 
requirement is that the thermal conduction parameter 1 * 

Also, for a plasma in which the electron and ion temperatures 
are uncoupled, the equilibration parameter X « 1 and, if 
the electron temperature is large compared with the ion tem- 
perature, then >-6 e /e i « 1 is required. 

For the limit e , X \ X6/9. , •* 0 , the governing 

equations simplify considerably. The continuity equation (2-6) 
becomes &C/5 t ~ 0 so that the density profile remains 



30 

stationary C = C (y) • The. momentum equation (2-7) becomes 
6U/&T =o so that the velocity remains zero throughout the 
process. The electron entropy equation (2-8) becomes 



i 



(2-16) 


The radiative transfer equation (2-10) remains the same, as 
does the absorption coefficient (2-11) and all three equations 
of state (2-12, 13, 14) . The electron entropy equation (2-16) 
plus the radiative transfer equation (2-10) , the absorption 
coefficient (2-11) , and one equation of state (2-13) form a 
complete set of four equations in four unknowns. 

These can be reduced to a single equation for electron 
temperature. Solving (2-13) for , taking the time deriva- 
tive and eliminating 5cs /b T using (2-16) yields 

© 


b8e oK . 
br S 


(2-17) 


Then solving 


i « 


2 fl-Cl 1/2 
5 <x£ £>t 


5/2 


(2-18) 


Bifferentiating with respect to y 


1/2 bV /2 

bi ^ 2 (1-0 V e 


1 - i c 


b9 5y/2 
d£ e 




a£ 


6ybT " 5 aS 2 (i_ £ ) 1 / 2 &T 


Combining this equation with (2-18) using (2-10) to eliminate 
i and (2-11) to eliminate K yields 



(2-19) 



t> 2 e 5 / 2 

e 

6y t>T 


5 C(l-C) a y 
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be 5 / 2 


e 

bT 


+ 




which can immediately be integrated once on T to give 


>,q 5/2 -i if 2 

2 °e _ 2 1 2^ dC Q 5/2 + aC 

5 by 5 £(1-0 e (1-C) 1//2 


f(y) 


(2-20) 


f(y) arises out of integrating with respect to t and is 
determined by the initial conditions. If the initial condi- 
tion is B (y , 0) = B*(y) , then (2-20) becomes 
6 . ® 


o = 


2 _b_ 
5 by 


C e; 


5/2 _ e o5/2 


1 2 ~ 2^ dC r fi 5/2 

J “ 5 CU-O ^y L e 


®Stj c b . - e = 3 


gO 5/2 -j 
e 


(2-21) 


< 1-0 


This is the equation governing electron temperature for a 
stafcionary-frozen-nonconducting plasma. 

It is noteworthy that although the dependent variable 0 g 
is a function of both y , and t , (2-21) is essentially an 
ordinary differential equation in 8 g in that only the partial 
derivative with respect to y appears. The only precaution 
in considering (2-21) as an ordinary differential equation is 
that in the solution, arbitrary functions of t arise instead 
of constants (as in the case of a real ordinary differential 
equation) . 

In order to formulate a well posed problem, appropriate 
boundary and initial conditions must accompany the governing 
equation (2-21) . The initial condition is 
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■ 8>,0) = 8°(y) 
0 - 


(2-22) . 


The density profile C (y) remains the same throughout the 
stationary process, and hence is the same as the initial den- 
sity profile. 

The basic boundary condition is that the intensity at the 
outer edge of the plasma equals the intensity delivered by the 
laser, 

i (0 ,t) = df/d-r 


But a boundary condition on temperature 0 e must be developed. 
(2-18) is rewritten as 


be 


5/2 

e 

bT 


5. UC 
2 


i 


and the boundary condition is 


oe 


5/2 | 

&T 


_ 5 aC (0) di 

lv=0 2 Jl - C{0) , dT 


(2-23) 


If C <y) is such that C (0) = 0 , then the appropriate boun- 
dary condition is on the first derivative with respect to y 



_e 

bybT 


y=° 


5 dC d£ 

2 dy y=0 d-r 


(2-24) 


for C CO) = 0 , and a£/dy| y=0 / 0 . It is necessary to resort 
to the derivative in this case because (2-23) will turn out to 
be satisfied automatically when C CO) — 0 . The same basic 
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boundary condition on the intensity governs both cases. Equa- 
tion (2-24) could also be found by combining (2-18) and (2-19) 
at the edge of the plasma. (2-22) and (2-23) or (2—24) together 
with the governing differential equation (2-21) form a well 
posed problem. 

2. Analytic Solution 

A study of the governing equations shows that the nature 
of the heating depends on the size of the thickness parameter 
a. . a large corresponds to ^n g sufficiently small such that 
the underdense portion of the plasma is optically thiclc. The 
result is a "bleaching wave" process whereby successive layers 
are heated until they become nearly transparent , allowing the 
beam to penetrate to a deeper layer, a small corresponds 
to Vn e sufficiently large such that the underdense portion 
of the plasma is optically thin. Then, the whole underdense 
region is heated simultaneously. If a is sufficiently small, 
the laser will have only a very small heating effect on the 
plasma and would be ineffective in the production of a thermo- 
nuclear plasma. 

If the initial temperature of the plasma is much less 

than T then another small parameter can be defined, 
o 

0° (y) = ug(y) (2-25) 

where .the maximum value of g(y) is chosen to be one and p 
is the small parameter, p is called the initial temperature 
parameter since it is a measure of the initial' temperature. 
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The governing equation is now 


_ l r 

' 2 b / . 5/2 5/2 5/2 N 2 1 “ 2 C dC f 5/2 5/2 5/2 

0 = 5 b7 l 6 e 5 / " 5 { (1 - Cl dy V Q e ^ g 




aC 2 


(0 e -ug) 


(2-26) 


If the final electron temperature is in the neighborhood of 

T then 9 will be of order one. Hence, during the hotter 
o ' e 

stages of heating, the terms in p will be relatively small. 
If only terms of order one are considered, (2-26) becomes 
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3/2 &6 e 
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by ~ 5 Cd - O d y e 
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0 O + 0(u) 


(2-27) 


This is designated the "hot equation” since it corresponds to 
the hotter stages of the heating. 

' In the initial stages of the heating, 8 e will be near 
0° which is of order p . For this case 6 must be 

G ' 

rescaled; 8 e = pB^ where 8^ is of order one. Then, 

(2-26) becomes 
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( v 5/2 -* 5/2 ) - 
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d£ 

dy 
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a C 2 

7^ 


( V g) 


(2-28) 


This is designated the "cold equation" since it corresponds to 
the colder stages of the heating. 
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' Solving the problem for 0° « 1 involves solving (2-27) 

and (2-28) separately and then appropriately matching the two 
solutions to the boundary conditions, initial conditions, and 
to each other. Solving such a problem where different parts 
of the solution must be matched with each other requires the 
method of matched asymptotic expansions such as described in 
detail by Cole 

a. Solution to the Cold Equation ’ 

3 /2 

A significant parameter a/p appears in the cold 

equation (2-28) which may be large or small. There are two 

•a /2 3 /2 

Important cases, a » p. or a « , a thicker y or 

a thinner plasma, respectively. 

- For a thinner plasma a « p 3 ^ 2 , (2-28) becomes 


°-!£(v 5/2 - 5/2 )- 


§< 


cu - C) 


d£ 

dy 


(V 


5/2 


-g 5/2 > 


+ 0{a/p 3/2 ) 


Rearranging this equation gives 

£ *>. (v 5/2 -* 5/2 ) ■ £ w, ^4= . .(«/u 3/2 ) 

which is integrated 

0 « = r g 5 / 2 + -ilX^1 2/5 + 0(a/p 3/2 ) (2-29) 

£(t) is an arbitrary function of T that arises in the 
integration with respect to y . A( T ) is found by applying 
the boundary condition, which is (2-23) or (2-24). The result is 



which is integrable 


1 = ~pl <f 5 + a) 

where a is the constant of integration, a can be evaluated 
by applying the initial condition (2-22). Since $ (0) = 0 , 
a = 0 also, and with (2-25) , (2-29) becomes 

0(i * ( g 5/2 * f ^2 ) 5/2+ ota/li3/Z) ■ (2 - 30) 

the solution to the cold equation when a « \? /Z • Physically, 
this corresponds to the case where the entire underdense region 
of the plasma is simultaneously heated from the very beginning 
of the laser pulse. 

Foe a thicker plasma a » , (2-28) becomes 

0 = 2 JL ( 0 5/2 5/2) + (0 g) 


5 a by '-p a -IT 

2 1 ~ ^ (0 5 / 2 -g 5/2 ) 


(2-31) 


Since the derivative term is multiplied by the small parameter 
H 3 / 2 / a , the method of singular perturbations must be used for 
the solution. Thus, the solution to (2-31) is composed of an 
outer solution, © Q , where terms of order \i /a can be 
neglected, and an inner solution, ©.j^ , where the temperature 


changes rapidly in a very short distance. 
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The' outer solution is to order one, 

® 0 = g + 0(p 3/2 / tt ) . (2-32) 

Physically this represents the unheated parts of the thick 
plasma . 

Finding the inner solution requires expanding the scale 
of y , 


_ y- y 0 

y = - 37 2~ ; 

(t /a 


■where y D — y 0 ( T ) gives the location of the "transition layer". 
The other functions in (2-31), g(y) and C (y) , will he constant 
to order one in the scale of y . For example , 


u3/2 _ 

£(y) « C(y o + -~£- y) = c(y 0 ) 
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+ 0(fJ 3 /a- 2 ) 


Then to order one, the inner equation is 
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which can be integrated 
@ 3/2 d0 i 

' I '6 1 g(y ) = w(T) 


c (yj 3/2 

y + 0(u / /a) . 


Vi - G(y„ 


(2-33) 


m(T) ia the function of time that arises in the integration 
with respect to y 
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The integral on the left side of (2—33) can be evaluated 
using ordinary methods of integration. Conducting the integra- 
tion, (2-33) becomes 


3/2 


0 . 


1/2 


3 tf + 2 V } + 109 


1/2 

- 1 

il 1/2 “ 
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/a) . (2-34) 


The complicated expression on the left side (2-34) can be 
expanded for limiting values of 0^/g * 

For B^/s ' near one, (2-34) can be solved as an expansion 

for 0./g , 


0^ = g{l + 4 exp 


r e 


[g(y)3 3/2L 3 


^ + m(-r) - 


*/i - £ (y 0 ) 


y ■ + - - - 3 


+ 0(u 3/2 /a) 


which is valid for ®j_/9 " 1 i • 

fts y goes to positive infinity, 0^ is seen to reduce 
exponentially fast to g which is just the outer solution (2-32) . 
Thus, the inner solution (2-34) correctly matches the outer 
solution on the positive y side of the transition layer. 

For ®./g much larger than one, (2-34) can again be 
solved as an expansion, 

-2, . 2/3 

© = r 4 m(r) - \ y ° ■ — y 1 + 0(p 3//2 /a) , (2-35) 

1 L 2 2 Vi - c(y 0 ) J 
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which is valid for » 1 . © ± is seen to grow very large 

if y becomes large negative. Hence, for large negative y , 
(2-35) must be made to match the solution to the hot equation. 
This inner solution corresponds physically to the front of a 
bleaching wave where the temperature begins to rise signifi- 
cantly above its initial value. 

b. Solution to the Hot Equation 
When the energy in the laser pulse is large compared with 
the initial internal energy of the plasma, then generally 

(Vfinal 76 ® >>X * this case (V final ^essentially 
independent of 6°(y) . The "hot equation" (2-27) governs this 
case. Dividing (2-27) by gives a linear first order equa- 

tion in 8 e , 
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b is the constant of integration and will be determined by the 
matching with the cold solution. Hence, 
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2/52/5 - 5 
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(2-36) 


The expression (2-36) is only physically realistic when the term 
in brackets is greater than or equal to zero. It is seen that 
given a t , this term is positive for y less than a certain 

value. 

The physical impossibility of its being negative is seen 
if the resulting expression for i(y, 7 ) is written using (2-18), 
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The term in brackets becoming negative corresponds to the inten- 
sity becoming negative which is impossible. The hot solution 
is only valid for positive values of the term in brackets. The 
cold solution must be used elsewhere. Physically, (2-36) des- 

I 

cribes the temperatures behind the bleaching wave. 

*} 

if * 2/5 « 1 , (36) can be expanded with a / as the 
small parameter, 


aVV/5 
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+ 0(a 
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(2-37) 


This limit corresponds to the physical case of thin plasmas 
when the entire underdense region is heated simultaneously. 
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c. Hatching Hot and Cold Solutions 
In the various expressions that compose the solution 
(2-30), (2-32), (2-34), and (2-36), various ’.unknowns arose: 
b, m Ct) . and y (t) . These unknowns will be determined by 
matching the different components with each other for ‘the 
possible ranges of the parameter a . 


This limit corresponds to a very thin plasma. The cold 
solution is given by (2-305 . The first term in the paren- 
theses is of order one and the second term is of order 
S/2 

a/|i ' and is relatively small. Then, (2-30) can be expanded 
about g(y) (where g is non-vanishing) giving 


0 M " 


5 + 


a 
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c± 

1372 “ 




+ o(a/u 3/2 5 


• ( 2 - 36 ) 


It is immediately seen that at the end of the heating (5—1) , 

5 /2 

the term in a/p ' is still small compared with g and thus 
there is no need to seek a hot regime. The temperature 
remains in the cold regime throughout the laser pulse and is 
governed by (2-38) 

Physically, this means that very little of the laser energy 
is actually absorbed by the plasma. The whole underdense region 
is penetrated by the beam so that a reflected light wave will 
arise. Since the intensity of the reflected light will be imper- 
ceptibly diminished due to weak absorption, its effect will be 
to roughly double the heating due to the incide'nt light, i.e., 
the term in a/p 5 ^ 2 


is doubled. 
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Neglecting the effect of the reflected wave, the tempera- 
ture produced is 


6 

e 


® (y) + 




Ce ” (y ) ] 3/2 


+ 0{a/ M 1/2 ) 


This is shown schematically in Figure 6 for a linear density 

5 /2 

gradient, a constant initial temperature, and for a/\i — 1/10* 


This limit corresponds to a plasma sufficiently thich 
that significant heating occurs, but thin enough that the 
heating is done simultaneously even at the beginning. 

The cold solution is again given by (2-30) , but now the 
term in a/p'^'^ 2 dominates the term of order one. Expanding 
about the term in a (where § , £ are non-vanishing) gives 


2/5 C 2/5 (f i) 


2/5 


(1-0 


I7T 


+ 0 ( 1 ) 


(2-39) 


This must match with the small a limit of the hot solution 
(2-37) . The matching requires that the constant b in (2-37) 
be zero. Once again the entire underdense region is heated 
and a reflected light wave will arise whose intensity has very 
little attenuation since a is small. The effect of the 
reflected wave then is to double the temperature in (2—39) . 

■ Neglecting the effect of the reflected wave, a composite 
solution can be constructed which possesses the features of 
the hot solution (2-37) and cold solution (2-30) , 
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■ e = (e o5/2 + f -~=: + o(a/p 1/2 ) 

e v e 2 

\i 3/f2 « a 

This limit corresponds to a plasma that is heated at the 
beginning by a bleaching wave. The cold solution is given by 
an outer solution (2-32) and a transition layer solution (2-34) . 
As was seen before, the transition layer (located at y o ( t) ) 
was matched with the outer solution for y > y 0 (t) * Now (2-34) 
must be matched with the hot solution (2-36) . The asymptotic 
limit for ©^ large is given by (2-35) . The limit of 0 Q 
small in the hot solution must be calculated. This matching 
problem is demonstrated schematically in Figure 7. 9 e in 
(2-36) achieves small values near the point where the term m 
brackets vanishes. Given t , the term in brackets vanishes 
at some y w = y w (r) . Expand the scale of y - in the neighbor- 
hood of y w , 


y-y, 


= 




where v is a small parameter to be determined in the match 
ing. Then (2-36) becomes 


r 3 

= 1 4ov 

u 2 </i-c(y w > 


2/3 

(~y v ) J + o(p) 


This can be matched to (2-35) if v = p 3/2 /a , m(T) =0 , and 

y (t) = y w (t) , in which case y v = Y • 

A composite solution valid over all y cannot be written 

in this case but the components can be summarized, 



45 


5/a 

0 L»)l 



Figure ?, Hatching of the hot solution 
transition layer and the cold solution 
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a 2 'V /5 r,5 3/5 


0 - S75 I <f 5) 

e (l-C ) 1/5 L 2 


2/5 r y C 7/5 ay “| 2/3 

; s ° J ■ 


0(M) 


(2-40) 


for y < y w (T) , 


8 53 6 ~° + 0(p 5//2 /a) for y > y (t) , 

G 6 W 


(2-41) 


r „ -I 3 / 2 hr h n 3/2 «- 8 e ^ 1/2 
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= 0 (p 3/2 /a) 


- a C 2 (y J 




= ty-y w 3 + o(u /a) 


(2-42) 


A bleaching wave is seen to arise which separates hot parts 
of the plasma from cold parts. The location of the wave. 


y w (T) is given implicitly by 


3/S 


5) ' '« | a 2/5 J w 


y„(T) c7/5, 


C' /D av 


(l-C) 


1/5 


(2-43) 


As 5 increases from zero, y is seen to increase monoton- 
ically with * . The velocity of the bleaching wave can be 
found by taking the implicit derivative of (2-43) , 


fls _ -2/5 1/5 / 7/5 

v w= §<!“*> CX-C(y w )] /c(y w ) 

1/5 

For a, laser pulse with i ~ r , V ~ r for a constant den- 
sity. A positive density gradient has the effect of limiting 
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the maximum velocity of the wave and eventually causing it to 
decrease. 

If c. is very small, the wave speed approaches the speed 
of light and reaches the critical density very early in the 
laser pulse. From that time on, the laser heats the entire 
underdense portion of the plasma simultaneously and. the tem- 
perature profile produced is 

2/5 

« e *(f^|) + Ofa 475 ^) (2-44) 

h "hot spot" is generated near Q = 1 , the critical density. 

Of course, the original expression for absorption coefficient 
(2-11) breaks down very near £ = 1 and .remains finite at the 
critical density. Nevertheless, a hot spot is produced in this 
region. This limit is called the "simultaneous heating case' 1 . 
In this case the actual temperature produced will be roughly 
twice that given in (2-36) since the light wave is reflected 
at the critical density and travels backwards through the 
region with very little intensity attenuation (see 2-10} ) . 
Hence, most of the laser radiation is lost to reflection. 

If the thickness parameter a is larger than a certain 
value a = the wave travels much more slowly and does not even 
reach the critical density by the end of the laser pulse. a„ 
is a critical value of a that depends on the density profile, 
but is always of order one. Hence, in this case there is no 
isolated hot spot as for a small. This limit can be termed 
the "bleaching wave case" since a bleaching wave continues to 
move into the plasma throughout the laser pulse . The heating 
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is by no means “simultaneous". In this case all of the laser 
energy is absorbed. 

In summary, stationary- frozen-nonconducting heating of a 

cool plasma will occur in one of four regimes depending. on the 

sizes of the thickness parameter a and the initial temperature 

parameter p . If a « p 5/2 , then the plasma is only slightly 

heated and the heating is simultaneous. A weak hot spot occurs 

at the critical density and nearly all the laser energy is lost 

to reflection. If p 5/2 « a « p 3/2 , the plasma is heated to 

temperatures much hotter than the initial temperature but still 

roost of the laser energy is lost to reflection. The heating is 

3/2 

simultaneous and a stronger hot spot appears. If p « cl 
< a , a very rapid bleaching wave traverses the underdense 
region early in the pulse. After that the heating is simulta- 
neous and a significant hot spot appears at the critical density. 
During the simultaneous heating stage, a significant part of 
the' laser energy is lost to reflection but not as much as in 
the previous cases. If a > a c , a much slower bleaching wave 
••propagates into the plasma and the wave does not penetrate all 
the way to the critical density. No hot spot or reflected 
wave appears, and all of the laser energy is absorbed. 

' 3. Application to a Linear Density Profile 
The solutions given in the last section can easily be 
calculated for a linear density profile, C - Y • These results 
for three different values of the parameter a are shorn in 
Figures 8, 9, 10. For a linear density profile, a c = 5.2 . 
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Figure 8. Temperature distribution for a 
linear density gradient with « = 1000 
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nearly simultaneous heating 

electron 



Figure 9. Temperature distribution for a 
linear density profile with o<= 4 
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electron hot spo 



distance y 


Figure 10. Temperature distribution for a 
linear density profile with \x 5lz « o< « l 


Figure 8 clearly demonstrates the bleaching wave moving 
into the plasma. In this case (a = 1000} _ corresponding to an 
optically thick plasma, the bleaching wave only propagates 
about halfway across the underdense portion of the plasma by 
the end of the laser pulse. In Figure 9 (a = 4} corresponding 
to a medium thickness plasma, the bleaching wave reaches the 
critical density just before the end of the laser pulse, and 
a weak hot spot is generated at the critical density. In 
Figure 10 (a « 1) corresponding to an optically thin plasma, 
the bleaching wave travels almost infinitely fast in traver- 
sing the underdense region and hence is not shown. A strong 
hot spot arises at the critical density. 

While Figures 6, '9, ana 10 give the analytical features 
of the solutions, they are given in terms of dimensionless 
quantities and do not show clearly the temperatures and other 
conditions that are attained. Figures 11 and 12 are summary 
plots for laser wave lengths X^ = l-06p and 10. 6p , respec- 
tively. The “maximum temperatures’* given are the maximum 
temperatures without the hot spot. Naturally , when a hot 
spot appears, there will be a much hotter temperature in a 
very localized region. In the simultaneous heating part, 
the temperature is of order a 2/5 T 0 . In the bleaching wave 
part, it is of order a 1/3 T 0 since smaller densities are . 
heated. The expressions for the maximum temperature are, 

2/5 

T = 1.08 (J A) eV for y. « a < 5 .2 

o 
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T = 1.64 
e max 


x lO - ^ (J A) 1/2 (^n ) 1/6 X 0 V3 for a > 5.2 


(2-45) 


The "plasma length" is defined as the thickness of the 
plasma that is irradiated by the laser. This thickness will 
be R for the simultaneous heating case, but will be of order 
a ~l/6 R for the bleaching wave case, 


I. = R = 1.12 x 10 27 {7nJ~ 1 X 0 " 2 for P « a < 5-2 

(2-46) 

3.18 x 10 12 (J /A) 1/4 (Vn e )“ 7/12 X“ X/3 for a > 5.2 ; 

' CL 

where I» is in meters- 

The maximum density heated is equal to the dimensionless 
length heated for the case of a linear density gradient. Then 
the maximum density heated is one for the simultaneous part, 
and'1.53a' 1>/6 . for the bleaching wave part. These are not 

shown on the summary figures . 

It is noted that the 10. 6p radiation heats in the simul- 
taneous heating regime for much smaller electron density 
gradients than does the 1.06p radiation, i.e., a relatively 
long plasma may be in the simultaneous heating regime under 
10. 6p radiation but in the bleaching wave regime under 1.06p 
radiation. By comparing Figures 11 and 12, it is also seen 
that the temperatures achieved are independent of the laser 
wave length in the simultaneous heating regime. However, the 
longer wave length heating produces higher temperatures in 
the bleaching wave regime. 
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It has been previously noted that not all the radiation 
is absorbed in the simultaneous heating case but that much is 
lost 'to reflection. It is important to Know what fraction of 
energy is absorbed— both from the standpoint of efficiency and 
from the fact that reflected radiation may damage the laser, 
integrating the temperature {2-44} for the thickness parameter 
a « 1 over the width of the underdense region, the fraction 
of energy absorbed can be calculated. This fraction is doubled 
to take into account the absorption of reflected light, and 
the result is 



2(f a) 


2/5 


£ 




which for a linear density profile reduces to 


= 1;72 a 2 / 5 
J o 

These results are shown in Figure 13. For a > 5.2- there is 
total absorption and for o « 1 the absorption is given by 
the above equation. Intermediate values of a are represented 
by the dashed line which is an interpolation. 

It has been shown recently that another effect may nearly 
eliminate. the reflected light for sufficiently high intensities 
This is the anomalous absorption effect that occurs near the 
critical density. This effect was studied theoretically by 
Kruer et.al., 16 and is confirmed by an earlier experiment by 

17 

Golcker and Sizukhin. 
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Figure 13. Fraction of laser energy absorbed 
for a linear density profile plasna 


. se 

4. Summary 

Analytic solutions have been developed for the laser heat- . 
ing of a one-dimensional plasma under the assumption that the 
plasma is stationary, non-thermally conducting, and that only 
the electrons are heated (the ions being frozen at their 
initial temperature). .These assumptions were found to be 
reasonable for many cases of interest. In these cases, the 
laser pulse is sufficiently short so that the processes of 
plasma motion, thermal conduction, and electron-ion equili- 
bration do not have time to make a significant effect. These 
assumptions reduce the problem (aside from equations of state) 
to a pair of equations? a radiative transfer equation and an 
equation for the temperature containing an energy source terra.’ 
These equations were combined and solved for all ranges of * 
the thickness parameter a by the method of matched asymptotic 
expansions. 

The use of a one-dimensional plasma was motivated by the 
need to make the mathematics tractable. No real plasma is 
strictly one-dimensional but many plasma geometries of interest 
are nearly one-dimensional and the analytic results of this 
study can be applied to such plasmas with reasonable accuracy. 
An example of thermonuclear interest is the longitudinal heat- 
ing of a column of plasma that has a uniform density 'across 
the column, such as might be found in a plasma focus or a 
theta pinch. Another example is the radial heating of such a 
column which is quasi-one-dimensional if the absorption length 

1/S is less than the radius of the column. Such a plasma 

o 

with a number density of 10 19 /cm 3 heated by a longitudinally 
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aimed laser of wave length 10. 6u to a temperature of 10 KeV 
could be contained with a magnetic' field -of 3 MG (as might be 
developed in a plasma focus) . For plasmas that are not nearly 
one -dimensional, the results of this study are still valuable 
in understanding qualitative results. 

The analytic solution shows markedly different behavior 
for different ranges of the thickness parameter a . For 
thick plasmas (corresponding to a large) the response is 
characterized by a bleaching wave proceeding into the plasma. 
Behind the wave is a hot nearly transparent plasma being irrad- 
iated by the laser. Ahead of the wave is an unheated opaque 
plasma. Physically, the motion of the wave corresponds to 
the heating of successive layers of opaque plasma to tempera- 
tures at which they become transparent and allow the radiation 
to pass onto deeper layers. For a greater than a certain 

distinct value a which depends on the density profile, the 
' c 

bleaching wave fails to penetrate to the critical density by 
the end of the laser pulse. Then none of the radiation is 
reflected back at the critical density. 

For a less than a c , the bleaching wave reaches the 
critical density before the laser pulse ends. At that time, 
reflection occurs at the critical density and some of the light 
is reflected. Although the heating is simultaneous once the 
bleaching wave traverses the underdense region, it is by no 
means uniform — an extreme hot spot is created at the critical 
density point. For relatively thin plasmas (a somewhat less 
than one) the bleaching wave traverses the underdense region 
extremely fast and all but the earliest stages of the heating 


are "simultaneous”. 
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For the thickness- parameter ' a smaller than the small 
parameter f , the entire heating process is simultaneous 
Yfifch no bleaching wave at all* is the initial temperature 

parameter. The efficiency of the heating drops further as 
more' of the laser energy is lost to reflection. For a even 
smaller (less than ) , the heating is so ineffective that 

the temperature is scarcely changed at all. 

It is clear that in the preparation of the plasma to be 
heated, the temperatures (or other conditions sought) can be 
optimized by proper choice of the thickness parameter a . 

For a given laser energy per unit area and wave length, cl 
depends on the density gradient. Hence, proper tailoring of 
the density gradient beforehand will give the best results, 
whether that be maximum temperature or maximum yield from a 
thermonuclear reaction. 

The solution found in this section is only for a stationary, 
nonconducting, frozen ion plasma. But the solutions for a 
stationary- frozen-conducting plasma, and a stationary-equili- 
brium-nonconducting plasma are closely related. In the sta- 
tionary— frozen— conducting regime, the effect of conduction will 
simply be to smooth out the nonuniformities in temperature that 
.are produced. This case is studied in more detail in the next 
section. In the stationary-equilibrium-nonconducting regime, 
the effect of the equilibration will be to divide the energy 
equally between the electrons and ions. In this case, the 
result is almost identical to the frozen ion solution presented 
in this section. 
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C» STATIONARY -FROZEN-CONDUCTING PLASMAS 


For sufficiently thin plasmas, thermal -conduction effects 
become important. When thermal conduction is strong enough to 
smooth out major temperature nonuniformities , then one has the 
stationary-frozen-conducting regime. This regime is similar 
to that studied in the last section except that thermal con- 
duction effects are dominant. In this section, a simple solu- 
tion is developed for the stationary-frozen-conducting regime. 
It is not intended that this solution be as rigorous as that 
of section B since it is desired only to demonstrate its 
major features. 


1 . Equations 

The requirements for this regime are e « 1 , X « 1 , 

X6 /8. « 1 as before, but now a different condition » 1 
© i 

Ion conductivity ri^ I s always small compared to electron 
conductivity and is neglected. It is assumed that e , X , 
10^/01 are actually zero, and n e » 1 • Then as before, 
the solutions to the continuity and momentum equations are 
q = £(y) and U = 6 , respectively. The electron entropy 
equation (2-8) becomes 


Do 

e 

i>^ 


. ttK l 
“ £ 


+ T) 


*V /2 

2/7 ° °e 


e C8 


hy" 


This can be combined with an equation of state (2—13) to 
eliminate the electron entropy, 
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C 



= OKx + % r\ 
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sV / 2 
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Using the absorption coefficient (2-11) this becomes, 


c 



o^u-o 


lt , ‘ 2 »r 

m L 'l "e 


(2-47) 


In the stationary- frozen-conducting regime, the thickness 
parameter cl is always small (see figures 2-5) . The sta- 
tionary-frozen-nonconducting solution has the temperature 
of order cl 2 ^ 5 ' as seen in the last section (in the simulta- 
neous part of the regime) . Since the primary effect of thermal 
conduction will be to smooth temperature nonuniformities, it 
is reasonable that the same temperature scale would appear in 
the adjacent conducting regime. Taking 9 e = a 0 a , (2-47) 
becomes 


68 c 2 &2 0a /2 

c ITT = e 2/2 (1 _ 0 J^ 1 + 7 ai1 e “^2” 


(2-48) 


If an <<: 1 the conduction term is negligible and the solu- 
‘e 

tion is stationary-frozen-nonconducting* If ar\ Q >:> 1 then 
thermal conduction is dominant and the solution is stationary— 
frozen-conducting. Thus it is seen that cit^ (not n e ) is the 
important parameter that determines whether the plasma is 
thermally conducting or not. This will doubtless shift the 
boundary between these two regimes as will be discussed in 


section D. 



But in the 
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Two unknowns remain in {2-48} ; © a and i . 
stationary- frozen-conducting regime, a « 1 so the laser 
intensity is only weakly diminished in passing through the 
plasma. The intensity is i = $ ' (t) [ 1 +■ 0(a / )] as seen by 
comparing the radiative transfer equation (2-10) and the absorp- 
tion coefficient (2-11} . In view of this fact (2-48} becomes. 


6 2 6 y 2 68 

° «\L £ Ot 


C 4 § (t) 

%i /2 \i~a x/2 


•] - 


2/5^ 




(2-49) 


Still needed to form a well posed problem are the boundary 
and initial conditions. It is not the intent of this discus- 
sion to include the details of finite initial temperature as 
done in section B for the nonconducting solution, so the initial 
condition is taken to he 9 (y,0) » 0 . The boundary conditions 
are on the -thermal conduction at both ends of the plasma, since 
there is a vacuum at the front end, y = 0 , there can be no 

N, 

thermal conduction, so 60 /6y| =0 . The same condition 

would apply at the rear edge of the plasma. 


2. Analytic Solution 

Consider the asymptotic expansion, = + + 

Then (2-49) becomes 


2 fe20 J /2 ’ , S 5 / 2 0 7 60 o 

«le-L by 2 &T 
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The equation to order one is b 2 0^ /by = 0 and has solutions 

6 = [a x (T)y + a 2 (T)} 2/7 where a 1 and a 2 are the arbitrary 

functions of t that arise in the integration with respect to 
y . The adiabatic boundary condition at y = 0 requires that 
a ^( T ) = o so that 8 q becomes purely a function of t . The 

equation to order l/o.r\ Q then becomes. 


0 = 


e 5/2 




C 2 ? • (f) 

8 3/2 ( i: £ ;V2 


Integrating once and applying the boundary condition at y - 0 
to order l/ari e yields. 


0 


,5/2 
o by 



C 2 dy 
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Integrating again and solving for 0-^ yields, 


9 i - -^(^r0^ ay } ay ■ 


+ a 3 (T), (2-50) 


where a 3 (r) is the arbitrary function of time that arises in 


the integration with respect to y . Then, 


0 = 0 (t) + - 

a 0 ari 6 


“o • 


,2/5 
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a. Underdense plasmas 

First suppose that there is no thermal conduction through 
the point y = 1 . This is the case where the plasma doesn't 
achieve critical density anywhere and y — 1 is the rear 
edge. Then the rear boundary condition is 60 a /by | y= j= 0 * 
This is satisfied automatically to order one. The order 
1/a.ri equation gives , 
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Integrating with respect to time, and applying the initial con 
dition of aero temperature. 


-Lf • J 


c 2 dy 

o (i-C ) 1/2 LJ o 




,-n 2 /s 


(2-51) 


Thus it is seen that applying the rear boundary condition 
to order l/ari e determines the unknown function 0 q (t) in 
the term to order one. Applying the rear boundary condition 
to the next higher order would determine the unknown function 
a 3 (T) in the term to order l/ari e . 

b. Overdense Plasmas 

Suppose that the plasma does achieve critical density and 
that the portion beyond the critical density (the dark part) 
is heated by thermal conduction from the irradiated portion 
■which is heated by absorption. This is shown schematically in 
figure 14. If the mass of the dark part is of the same order 
of magnitude as the mass of the light part, then the presence 
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o £ the dark part will not alter the hasic character of the 
solution, i.e., that it is dominated by thermal conduction. 
However, the temperature produced will t>s less since a 
greater mass of plasma is being heated. 

To construct the solution in this case, one must solve 


(2-46) for both the light part, 6^ , and dark part 8 ad , 

separately. Then the two solutions are matched at the point 
y = 1 with the condition that thermal conduction at the 
junction is the same for both. Thus {0 a £- e ad > l y= i = 0 and 
(b0 ajJ /&y - bS^/by) [ 0 . The dark form of (2-49) is 

the same as the light form except that the absorption term 
(that containing k ' (t) ) is missing in the dark equation. 


The order one dark equation is the same as that for the 

light equation; &^oJ d ^/hy^ = 0 , and its solution is 

e = [b (T)y + b (t)] 7 "^ 2 . Applying the adiabatic boundary 
od X 2 

condition at the rear edge (y = L) of the plasma, bO^/by}^- 

requires that b^ (r) =0 so that 0 Qd = ® od ( T ) - Matching 

this to the light solution to order one requires that 

0 , (t) = 0 „{t) =6 (t) . The order l/an equation is, 
od o£ o ® 
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Integrating with respect to y and applying the boundary 
condition at y = L gives, 
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Integrating again gives, 

1 W . 
6 ld - 77? TF 
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f-f{fW}ay + b 3 (T) , 


y y 


where b 3 (t) is the arbitrary function arising in the integra- 
tion. Matching to the light solution (2-50), ~ 

e ue (1 ’’ rI merely forms a relation between the two unknown 
functions of time a 3 (T> and b^r) . But the matching of 
the first derivatives &0 ld / b yl y;=:L = b0 i£ /by ly=l glVeS> 
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integrating with respect to t and applying the zero initial 
temperature condition, 

»o - [ f 5 C :d^72 {£ T 


Which is the same as the solution to the previously solved 
case except for the different range of integration in the 
denominator (over the whole plasma) . 

Then the usual nondimensional electron temperature is, 
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(2-52) 
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Thus, for the stationary- frozen-conducting plasma, the 
temperature . is seen to rise — uniformly — with the two-fifths 
power of the laser energy released (as measured by . This 
agrees with the fact that thermal conduction dominates and 
also with the fact that as the plasma heats up, it becomes 
more and more transparent to the radiation, thus heating 
more and more slowly. The fall off in the rate of heating 
is clearly seen in figure 15 which is the time history of the 

electron temperature 0 e for a plasma with a parabolic den- 

3 12 . 

sity profile, £ = y - j y , heated by a constant intensity 

laser pulse for which « = t . This solution is based on the 

assumption that the masses of the light and dark parts of the 

plasma are of the same order. 

Suppose that the dark part of the plasma is much larger 
than the light part. Then, even though the light part is 
dominated by thermal conduction and has a nearly uniform tem- 
perature, the dark part may be quite nonuniform. If>the dark 
part is large enough, then it will have a nonuniform tempera- 
ture profile diffusing with time into the dark part. Since 
the thermal conduction coefficient of a fully ionized plasma 
is proportional to the five-halves power of the temperature, 
the thermal diffusion is nonlinear. It will proceed into the 
cold dark part as a "thermal wave" resembling the "bleaching 
wave" found in the stationary-frozen-nonconducting regime. 

The close relation of these two phenomena is discussed in 
detail in chapter IV. 



70 


electron 
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Figure 15. Time history of electron temperature of a 
plasma in the stationary-frozen-conducting regime 
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B„ RE-EXAMINATION OF HEATING REGIMES 

As was previously mentioned, it is necessary to take any 

solution found and reconsider the regimes in which it is valid. 

This is especially true in this study since the dimensionless 

electron temperature 9 e " varies over many orders of magnitude. 

6 is of order a’ / ' 5 in the stationary- frozen-conducting 
e 

regime and in the simultaneous heating part of the stationary- 
frozen-nonconducting regime. For the bleaching wave portion 

Q 1/3 

of the stationary-frozen-nonconductrng regime a ~ a 
a is a parameter that may either be much larger or much 
smaller than one. 

Any locally extreme values that may occur in 9 g or its 

derivatives may also invalidate the solutions. Indeed such a 

singularity does occur in the neighborhood of the critical 

density where a temperature singularity does occur (in the 

stationary- frozen-nonconducting regime). However, it is not 

the purpose of this discussion to probe the fine points of 

local phenomena — such as that very near the critical density. 

Rather it is intended to re-examine the overall behavior and 

to find the correct boundaries of the various regimes since 

0 is not of order one. 
e 

1. Onset of Motion 

The original estimate of the plasma regimes (section A 3) 
indicated that the onset of motion occurs toward the upper 
right corner of the summary charts (figures 2-5), i.e., for 
higher values of ?n e and J 0 /A . A study of the governing 
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equations with cognisance of the correct temperature scale 
reveals that the original estimate is too conservative, and 
that the onset of motion occurs for even higher Vn e . 

The appropriate criterion to determine the onset of motion 
is the conditions for which the plasma motion begins to sig- 
nificantly -affect the temperature that is produced. Consider 
the asymptotic expansions. 


C = v 0 < a >£ 0 + v x ta > G >?l + '•» 

U = ‘ ^(a, e)!^ + ... 

0 = 1 (ct) 0 + X 7 (<x, e) 0 + • « « (2—53) 

6 o o J- 

1 = i + 6 1 (a,e)i 1 + ... 

y = 0y 


where X (a) and v (a) are the orders of magnitude of the 
o o 

^temperature and density for the stationary solution. 0 is the 
appropriate length scale. They are to be determined by the 
nature of the stationary solution. The terms of order 


% \ h contain the lowest order effect of motion. 
The terms of order one , v Q and 9 ive the stationary 
solution. The object of this discussion is to find X^A 


which is the ratio of motion induced temperature change to the 


stationary temperature. 

Applying the expansions (2-53) to the momentum equation 
(2-7) together with the equation of state (2-12) gives 



?3 


[ 0 ] + [hw 1 + nr'j - ^| {c * 9o) ] v ** • = 0 * (2 ~ S4) 

In •(•his and in subsequent equations, the first bracketed 
term represents the stationary terms and the second represents 
the lowest order effect of motion. 

Applying the expansions (2-53) to the continuity equation 
{2-6) gives 


IS&7 1 ] + + 


6g i 


P 




4 • . . — 0 


(2-55) 


Applying the expansions (2-S3) to the radiative transfer 
equation (2-10) with the absorption coefficient (2-11) gives 


ri 

\J & y + y72 "P72 


I 2 
v o a 


Co 2 i 0 


f 




00 •/! “ v 0 S, 


1 , f 6 l. &1 l . 6 l v o a .?e i i . 

=1 LT 7~" 737? 772 + 


6y 


Aq 


V(j 2 v 0 ^ O i 0 8j -j 

”7372* n 3/2 “ 7 ^ 72 ' , fl 5/2 J 

A 0 0 C A 0 20 e 


+ — 0 


{2-56) 


Note that the radical Jl~~ \> Q ^ 0 +'J 'l' was not expanded. 'This 
factor becomes important only if the plasma becomes overdense 
at some point. 

* (The electron energy equation is found by combining (2-8) 
With (2-11) and (2-13) , 


C 


550 

e 

5t 


20 


■K 

Or 


£u(e 


!!s 

Qy 


20 


+ thermal 
conduction 


e . j>£\ _ 

3 0y/ “73/2 1 — ; and equili 

°e ' S ' A “'> b ration terms 
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The thermal conduction term is not written out since its only 
effect is to smooth out thermal nonuniformities . The equili- 
bration term is not written since it can only affect the 
electron temperature by a factor of two. Applying the ex- 
pansions (2-53) to the electron energy equation gives 


£ v oAo(,C 


°5t 3 9 0 5t 


&e< 


2 . 

Co x o 




=l 


. he 2 b^ 1 f 6s 0 2 

[voAi(c 0 &r + ) + v 1 a °UiSt 3 0 °&r-) + 


§l v oA 0 e 

p 




2 

6y 


\ , v " v l a 2 Cogl io . 
) “772* 0 V2 


Ao 


2 2 . 2 2 . 
VqB^CX. V 0 


Ao 


372-^72 


l l + _o (2-57) 

572 2 ^/2 } + **• " ° * (2 5?) 


The most general equations are found when all terms 
within a pair of brackets are the same order of magnitude. 
Then from the stationary terms, one can conclude that 


Ao 


3/2 


Pv 0 

. 5/2 


= a 


Vo 


a . 


(2-58) 


From the motion terms, one can conclude that 




75 


v 0 Ao 2 

~T 


h ‘ 


'o 2 


Al 


2 

Ao 2 

7 e - 


(2-59) 


(2-58) and ( 2 - 59 ) are now examined for various cases to 
determine the appropriate variables for the onset of motion. 
For a linear gradient plasma in the bleaching wave 

regime, v 0 = B = a 1/5 and A 0 = <* 1/3 (from section c > 
which is consistent with (2—58) . Then (2-59) leads to 


1/2 1/22 2/3 2 

gj = a e , vp = a e , 6 -j - a e , ana Al 

Thus Aj/Aq = a 2 // 3 e 2 is the parameter determining the onset 
of motion in this regime. 

For a linear gradient plasma in the simultaneous heat- 
ing regime, Vo - p = 1 and Ao = a 2/5 • This is consistent 
with the second equation in (2-58) but not the first. The 
reason is that the characteristic length scale predicted by 
the radiative transfer equation is the absorption length, 
which is erroneous in this case since the characteristic 
length is the distance to the critical density where the heat- 


o.e 


2/5 , 


ing end 3 . Then (2-59) leads to 5 . - a e 

4/5 2 

Al “ a e 


v 2/5 2 , 
V ! = a e 


, _ 2/5 2 

6 X - a e 


Then Aj/Ap = a 2 / 5 e 2 is the 


parameter determining the onset of motion in this regime. 

is the specified scale of the 


For a uniform plasma, v 0 
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"density (which is constant) . Then (2-58) leads to = 

(av 0 ) ^ // ^ . * P = a 2 ' /5 v 0 "* Thus (2-59) leads to = 

4/5 9/5 _ 6/5 21/5 2 . 6/5 16/5 2 

a v 0 , e , - a v 0 e » 6 1 " a v o e , and 

AJ. = a 8/5 Vo 18/3 e 2 . Then hj / ko = a 6/S v 0 16/5 £ 2 is the para- 
meter determining the onset of motion for a uniform plasma , 

In a linear gradient plasma, the onset of motion will 
occur in the simultaneous heating regime for cases of interest. 
In this regime the important parameter is 

6* = € 2 <x 2/S ■= 2 -42 x 10" 47 (J /A) 2/S (7n ) 2 X 4 t 2 . (2-60) 

O 6 ® P 

2 2/5 

If e <x « I , then there is no motion during the laser 
2 2/5 

pulse. If e a ' » 1 , significant motion occurs. 

In a uniform plasma, the important parameter is 

«* - a* /5 v“'V - 1.20 x 10 24 (d o A)- 4/5 c 16/5 C^ . < 2 - 61 > 


where £ has replaced V(} since they are identical for a 
uniform plasma. (2-61) fails very near q = 1 since the 
radical £ "'+ 77.' becomes very small in that case. 


2. Onset of Thermal Conduction 
In developing the stationary-frozen-conducfcing solution 
for a < 5,2 (section C) , it was found that aT , e is the 
parameter measuring the thermal conduction. However, thermal 
conduction will first become dominant in regions where a > 5.2 
and a different parameter is required to determine the onset 
of thermal conduction. 

The energy equation including thermal conduction but not 
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-the motion or equilibration terms, is 


be. 


c 2 i 


.2,7/2 


tT8 


V /2 /I^ + 7^ &y 


This was found by combining {2-8} , (2-11) , and (2-13) . With 
the scaling C ~ v 0 Co > 6 e = A ° 9 ° ' y = ^ > the ener 9Y 
equation becomes _ 


be. 


v v oa„ 


Coi 


" v oCo 


7 


5 /2 

The scales are consistent if the parameter v 0 a / A 0 = 1 • 

Then onset of thermal conduction occurs when the parameter 
Ao 5 / 2 rie /8 2 Vo exceeds one. Eliminating Vq , the ke Y Para- 
meter is ari e /p 2 * 

For a linear gradient plasma in the bleaching wave regime 
£ = ^ a _1 / 6 and the onset of thermal conduction is determined 
by the parameter - 

a 4 / 3 ^ = 3.07 x 10" 30 (a o /A) 1/2 (7n e ) 7/6 X o 8/3 t p . (2-62) 

If a 4/3 \ « l , there is no thermal conduction, and if 
g. 4 / 3 ^ » 1 , thermal conduction is dominant. 

For a uniform plasma, consideration of (2-40) indicates 
that the plasma length heated is Y max ~ 1*13 a 2//5 C > 

valid for q less than one. Then p = cT 2 ^ 7//5 is the 
appropriate choice and the thermal conduction parameter is 

J « a V5 p e C 14/S =■ 5.51 x 10 11 (J o A)- 1/S C 14/S X“ 2 t p . (2-63) 
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3. Onset of Electron-Ion Equilibration 
The electron and ion energy equations including equilibra- 
tion but not thermal conduction and motion terras, are 




2 . 2 

n •fVj' "*- — X 3'72 ^ ** ®i^ * 


+ X T372 {e e ” e i ) 
e 


(2-64) 


With the scaling q = v 0 £o , 9 = Ao e o , = Ao 9e>£ , y - Py 

3 /2 

the scale of the equilibration terms in (2-64) becomes vo^/Ao 
which is the proper scaling when the ion temperature is of the 
same order as the electron temperature. If however, the ions 
are very cold (at the initial temperature) then u is the 
proper scale for 0 i . In this case the coefficient of the 
equilibration term in the ion energy equation becomes - V(1 1 / a/ U 


which determines the ion heating. 

* 3/2 

.-The parameter X eq = v 0 X/a d determines the bounds of 

* 

the equilibrium regime. If » 1 there is electron-ron 

* 

equilibrium. If X « 1 the electron and ion temperatures 

. * X/2 

are not in equilibrium.- The parameter X fr = v 0 X/a 0 U 

* 

determines the bounds of the frozen regime. If X fr » 1 
some ion heating takes place, though equilibrium is not assured. 
If « 1 the ions remain frozen at their initial temper- 

ature throughout the heating. 

Before calculating these parameters it is necessary to 
know the expression for the small temperature parameter 
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H = T°/T q . Using (2-1) , 

W = 3.0 x 10 35 (J o /A)" 1 (7n e )~ 1 X o _4 T° , (2-65) 

where T° is in eV and the other terms are in the usual units. 

In a linear gradient plasma, in the bleaching wave regime 
1/3 

A 0 =* n so that the equilibration parameter is 

X eq “ X/a 


2/3 


7.05 x 10 7 (J A)~ 1//2 (vn ) 1/6 X~ 4/3 t 
o e o • p 


for a > 5.2 , (2-66) 


In the simultaneous heating regime where j\ = a 


2/5 


X* = X/a 3//5 = 5.48 x 10 13 (J /A) -3 / 5 x“ 2 t for a < 5.2 

eq ° ° P (2-67) 

The freezing parameter in the bleaching wave regime is 

Xf = X/ a 1/3 U = 67.2 (Vn) 1/3 X” 2/3 T° ** for a > 5.2 

tr e o e p (2-68) 

In the simultaneous heating regime the freezing parameter is 
. * 

X” r = X/ a 1/5 u => 4.08 x 10 13 {J yA)-V 5 x-2 T o 


o e p 

for a < 5.2 (2-69) 


In a uniform plasma v 0 = C ana X = ( a C) 
equilibration parameter is 


2/5 


and the 


X* = 5.50 x 10 13 (J /A)- 3 / 5 £ 2 / 5 x" 2 t , 

eq ' o v op' 


(2-70) 


and the freezing parameter is 


-1 


X fr = 4 * 09 x 10 13 (J /A)“ 1/5 £ 4/5 X~ 2 T° t - 


(2-71) 



so 


4* Summary 

The results of calculating the various parameters are 
shown in figures 16, 17, 18, 19 for a linear gradient plasma, 
and figures 20, 21, 22, 23 for a uniform plasma. The lines 
bounding the regimes are extended on into the conducting 
regime as dashed lines even though these parameters were 
calculated from the nonconducting point of view. 

The heating of a linear gradient plasma falls conveniently 
/into several regimes. For short plasmas (high Vn & ) the 
r heating is nonstationary. For slightly lower values of 
the heating is stationary and conducting. For even lower Vn & 
{longer plasmas) the heating becomes stationary and noncon- 
ducting. Equilibration effects depend primarily on the energy 
deposition and the pulse time rather than on . The 
solution of section B applies in the stationary- frozen- 
conducting regime, and the solution of section C applies to 
the parts of the stationary-frozen regime above the line 
a ~ 5.2 . 

These solutions can also be applied to the regions where 

some," or complete equilibration occurs. This fact is easily 

seen by observing the stationary energy equations for electrons 

'and ions, (2-64). For the frozen case, the electron energy 

equation governs the heating as studied in sections B and c. 

For equilibrium regimes, = 6^ and the two equations combine 

to give , ? 

4>e e _ ' c i ^ 
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Figure 1 6 ,. Re-examined heating regimes for a 
nanosecond pulse from a Kd* glass laser 
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~6 

w = 10 sec X o = 10. 6y. 



Figure 18. Re-examined heating regimes for a 
microsecond pulse from a C0 2 gas laser 
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t « 10” 6 sec \ « 



Figure 22. Keating regimes with a uniform plasma 
for a microsecond pulse frota a CO, gas laser 
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—9 

t = 10 see \ = 10. 6u 
P 0 



Figure 23. Heating regimes with a uniform plasma 
for a nanosecond pulse from a COg gas laser 
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This is the same as the previous electron energy equation 
except that a/2 has replaced a as the absorption parameter. 
Hence, the equilibrium solution is the same as the frozen 
solution except a is replaced by a/2 . if there is partial 
equilibration, the absorption parameter will be between a/2 
and a . A similar calculation including the thermal conduc- 
tion term in the electron energy equation shows that in an 
equilibrium regime, the thermal conducting parameter is reduced 
by one half. Thus the solutions found in sections B and C 
apply to essentially all the nonstationary regimes. 

The heating of a uniform plasma falls into the same basic 
regimes. For a uniform plasma, larger values of £ corres- 
pond to shorter plasmas and vice versa. The length of the 
plasma is just the length that is heated. One noticeable 
difference of the uniform plasma from the linear gradient 
plasma is that even for large energy deposition, the heating 
may still be in an equilibrium regime for the longer pulse 
lengths . 

As before, the solution found in section B applies in 
the stationary-nonconducting regimes whether equilibrium, frozen 
or somewhere "in between. The solution of section C does not 
apply in the conducting regimes because there is no critical 
density point to cut off the radiation and limit its penetration. 

No solution was developed for those conducting regimes 
where the critical density does not appear {such as in the 
uniform plasma) or where the critical density region may not 
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be irradiated (such as in the part of the stationary-conduc- 
ting regimes below the a = 5.2 line for the linear gradient 
plasma) . It is expected that a wave like heating will occur 
in these regimes also, and the reasoning behind this is 
explained following. The thermal conduction will tend to move 
as a wave. Behind the wave is' a hot transparent plasma, and 
ahead of the wave is a cold opaque plasma. Then the laser 
beam will tend to penetrate the transparent hot region and 
be absorbed near the front of the wave. Thus, even though 
thermal conduction is strong, it is expected that the speed 
of the “wave" will be governed primarily by the absorption 
of radiation in its front, i.e. as in the bleaching wave. 

Hence it is expected that the solution in this regime will be 
very similar to the stationary- frozen-nonconducting solution— 
both insofar as the temperatures achieved and the lengths 
heated. Of course, detailed phenomena of the two solutions 
may differ. Thus the extension of the lines for the onset of 
motion, freezing, and equilibration into the conducting regime 
in figures 17-23 is apparently justified. 
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E. POST-HEATING DYNAMICS 

One objective of laser plasma heating is to achieve gain - 
in a thermonuclear reaction, The general criterion for gain 
is governed by Lawson's criterion which is discussed in 
Appendix 2. Lawson’s criterion requires that the product of 

the number density of the reactants (deuterium and tritium 

14 3 

ions) and the time of confinement must exceed 10 sec/cm , 
and the temperature must exceed 10 KeY. A certain minimum 
confinement time is required simply to get enough energy out 
of the reaction. For a CO^ gas laser, the appropriate density 

would be on the order of the critical density, 10 /cm . This 

4 * * * 

calls for a confinement time of 10 usee , For a Nd glass 

21 3 

laser, the appropriate density would be 10 70 s . This calls 
for a confinement time of 100 nsec. Thus, the plasma must 
be contained for times significantly longer than the laser 
pulse that would heat it in a stationary regime. Hence, what 
happens to the plasma after the pulse is of critical importance. 

■ 1. Post-Heating Dynamics for a Linear Gradient 
For a thermonuclear reaction to proceed at all, one con- 
dition must be met. The electrons must equilibrate with the 
ions before the plasma motion reduces the electron temperature 
and the densities. This is critically important because the 
thermonuclear reaction rate is strongly temperature dependent. 
Thus it it essential to find two time scales, the time scale 
for eiectron-ion equilibration, t ' , and the time scale for 
significant plasma motion to arise, t . 


A third scale of 
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interest is the time scale for significant thermal conduction. 


An examination of these time scales indicates that the 
post-heating dynamics falls into two regimes for conditions 
of interest. One regime is characterized by frozen-conducting 
dynamics, and the other by equilibrium-conducting dynamics. 

a. Acoustic Time 

The acoustic time is the time for an acoustic wave to 
traverse the heated region. The thickness of the heated region 


y R , and the maximum electron temperature. 


(which 


? 1/2 

determines the acoustic speed using a = 6,36 x ICr T q m/sec 
for T e in eV) , are given in section B. Then the acoustic 
time is 


2,54 x 10 11 {vn e )“ 2/3 l“ 2/3 sec 


(2-72) 


This is for the bleaching wave part, $ > 5.2 since it is in 

this regime that the three effect^ motion, equilibration and 

conduction, will compete, 

b. Electron-Ion Equilibration Time 

The correct electron-ion equilibration time corresponding 

' to the maximum temperature and density in a linear gradient 

£ 

plasma is given by * t /X (in the bleaching wave regime) 

eq p eq 

& 2/2 

where X * \/a * is equilibration parameter calculated 
eq 

in section D. Then, 

t eq = 1.42 x l<f 8 (J q /A) 1/2 (GnJ “ 1 / 6 ) _ o 4 / 3 


sec 


(2-73) 
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c. Thermal Conduction Time 

The thermal conduction time is the time for significant 
thermal conduction to diffuse accoss the thickness of the 
laser heated region and is given by 


_ p max C v^max R ^ 
th " K e 


(2-74) 


where K is the electron thermal conductivity given by 
e 

Spitzer.^ 4 Applying the maximum temperature and density found 
in section B 3, 


'th 


= 3.26 x 10 29 {J o /A)" 1/2 (7n e }“ 7/6 \' 8/3 


The same result would be attained by taking t fch - 
where is the thermal conduction parameter 
section. D. 


sec . (2-75) 

t /, 4/3 . 

tp/Ca r) e ) 

found in 


d. Summary 

The results of these calculations are presented in figures 
24, and 25, for a Nd + glass laser and a C0 2 gas laser respec- 
tively. It is seen that post-heating behavior falls mainly 
into two regimes: the equilibrium-nonconducting regime where 

t « t « t . . and the frozen-conducting regime for which 
eq a th' 

t « t « t -In the former regime, nonconducting 
til 3 6CJ 

dynamics will occur and in the latter the dynamics will be iso- 
thermal, based on the isothermal speed of sound. 

In the equilibrium-nonconducting regime, after the pulse, 
the temperatures equilibrate, and the subsequent motion is of 
an equilibrium nonconducting plasma. Thermal conduction never 
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occurs. 

In the frozen-conducting regime, after the end of the 
heating, strong thermal conduction smooths the nonuniformities 
in the temperature. Then the isothermal motion begins. Very 
little electron-ion equilibration occurs before the motion 
cools the 'plasma. At some' point in the cooling, the equili- 
bration time gets short enough that the ions equilibrate 
with the electrons — but at a much lower temperature than the 

original electron temperature. 

Thus the issue is clear, for a thermonuclear reaction to 
occur during the confinement time of the plasma, the plasma 
must be in the equilibrium-nonconducting regime, otherwise the 
ions never get heated to the temperatures required for a ther- 
monuclear reaction. This criterion requires thicker plasmas 
(lower Vn } , and the requirement is more stringent for grea- 
ter laser energy (larger J 0 /A) • 


2. 'Post-Heating Dynamics for a Uniform Plasma 
In this geometry as well as the linear density gradient, 
the relative sizes of the equilibration and acoustic times are 
all important. These and the thermal conduction time from 
(2-74) can be easily calculated using the maximum temperature 
and density in a uniform laser heated plasma. The same results 
would be found by taking t eg = t p /^ e q > an<3 fc th - 
The time scales are. 


t 

a 


8.66 x 10" 13 (J 0 A) 2/5 c 


-8/5, 2 
o 


sec 
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1.82 X icf 14 (cr A) 3 / 5 £~ 2 / 5 x _ 2 sec f 
o ° 



1.82 x 10~ 12 (J o A) 1//5 £“ 14/5 X o 2 sec , 

(2-76) 


The relative sizes of these time scales are presented in 

figure 26, along with the maximum temperature achieved (which 

is found by* applying (2-40) to a uniform plasma) . Also shown 

is the length of plasma heated, L corresponding to a lO.Sp. 

* + 

COj gas laser. The lengths for a l.OGp Nd glass laser are 

smaller by a factor of 100. The chart uses and 

£/(!-£) '^ as the independent variables. A nonlinear scale 

for C is shown side by side with that for £/{l-£) ' l As 

with the linear density gradient, the post-heating behavior 

falls basically into two regimes: the same two in fact. Given 

a constant J q /A , higher densities lead to frozen-conducting 

motion, and lower densities lead to equilibrium-nonconducting 

motion. This is the opposite of the natural tendency for t 

eg 

to decrease with increasing density. The reason for this 

behavior is that for higher densities, the temperature is 

3/2 

higher which tends to increase t^ (t ~ T e / n e ) * The 
temperature increase is greater than the density increase, 
hence the equilibriation time will be greater. 

a. Wave Diagrams for a Uniform Plasma 
Xt is appropriate to calculate a wave diagram for the 
motion that occurs in a uniform plasma. This calculation has 
been performed for a plasma in the equilibrium-nonconducting 
regime. This is the regime of interest for achieving a fusion 



froEen-conducting dynamics 



Figure 26 , Fost-heating dynamics for a uniform plasma 
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reaction, so this calculation will be applicable. To simplify 
the calculation, several assumptions were made: 1} neither 

thermonuclear energy release not bremsstrahlung , alters the 
dynamics, 2) the temperature profile produced by laser heating 
was taken to be uniform at the average heated value, T . The 
initial temperature was taken to b.e 1/100 'of T . The results 
of this calculation are presented in figure 27 which is the 
wave diagram, and figure 28 which shows the temperature 
profiles for successive times. 

It is seen that the plasma expands freely into the vacuum 
at the plasma edge. A shock tube like phenomena occurs at 
the inner edge of the heated region. An expansion wave moves 
into the hot plasma and a shock wave moves into the unheated . 
plasma. The shock wave heats the cold plasma to only 12% of 
the temperature of the laser heated region. 

A similar calculation was performed for a uniform plasma 
with a solid wall at the inner edge of the heated region. 
Physically this corresponds to a plasma heated from both ends, 
designed so that the bleaching waves just meet in the center 
at the end of the laser pulse. The wave diagram for this case 
is shown in figure 29. 

b. Release of Thermonuclear Energy 

The usual estimate of thermonuclear energy release is to 
assume the plasma reacts at its initial temperature for exactly 
one acoustic time. Of course this is not the case because 
the temperature begins to reduce from the very beginning, and 




Figure 2?. Ware diagram for a laser heated uniform plasma 
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Figure 28. Temperature profiles 


.8 


-.8 -.k o A 

distance y , b heated - re S ion 

plasma 
edge 

Figure 29. Ware diagram for a laser heated uniform plasma 
confined at the inner edge of the heated region 





Figure 30. Comparison of thermonuclear 
energy release to the usual estimate 


104 


it does not cool abruptly at exactly one acoustic time. It 
is worthwhile then to use the solutions of motion in figures 
27 and 29 to compare the actual release of thermonuclear energy 
to the estimate. 

This calculation was performed assuming the results of 
the wave diagrams. The result is shewn in figure 30. The ratio 
of the actual energy released by the D-T reaction per unit cross 
sectional area, E DT , to the estimate of the same, (E DT ) Q , 
is plotted versus the mean temperature to which the plasma is 
heated, T . The actual energy release is lower than the 
estimate by only by 50% at 10 KeV (or 20% for the case of a 
confining wall at the inner edge of the heated region) . ' 

4. Summary 

An examination has been conducted of the plasma dynamics 
following -laser heating in a stationary regime. Post heating 
dynamics will appear in one of two regimes. Frozen-conducting 
dynamics will occur for denser plasmas (since higher densities 
lead to higher temperatures) . In this regime, the electrons 
never equilibrate with the ions before motion dissipates and 
cools the plasma. No thermonuclear energy is released in this 
regime. For lower plasma densities, equilibrium-nonconducting 
dynamics occur. Here, the electrons equilibrate with the ions 
before the motion and thermonuclear energy is released if the 
initial temperature is high enough. The- usual estimates of 
thermonuclear energy released are accurate to within a factor 
of two or better for 10 KeV temperatures. 
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F. APPENDICES 
Appendix 1 

Effect of Electromagnetic Forces 
The primary thrust of this chapter is laser heating of a 
stationary plasma. Hence the effect of forces on the plasma 
does not appear in the solutions. Nevertheless, the size of 
the forces will affect the regime of validity of a stationary 
solution. The object of this section is to consider the size 
of the electromagnetic force that was included in neither the 
original equations, nor in any discussion of the regimes of 
validity. 

Electromagnetic forces are proportional to the intensity 
of radiation. Thermal forces (pressure) are proportional to 
the time integral of the intensity since absorption builds the 
temperature and pressure with time. Thus it is possible with 
a sufficiently high intensity to have an electromagnetic force 
larger than the thermal force. But even a large force may 
have little effect if it operates for only a very short time. 
The largest electromagnetic forces will arise in picosecond 
laser pulses — which have the shortest duration as well. Thus 
the ratio of the two forces does not determine their relative 
importance. One must find the -actual effect on the plasma 
motion. Previous investigations of electromagnetic forces have 
been made. Hora^’^ discussed the electromagnetic forces 
associated with the radiation and drev some conclusions relat- 
ing to its effect on' plasma motion. Lindl and Kaw^° studied 
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the effects of these forces when the radiation is not neces- 
sarily normal to the surface of the plasma. However, these 
two papers do not determine the relative effects of electro- 
magnetic and thermal forces for typical laser produced plasmas. 
In this appendix, it is shown that electromagnetic forces have 
little effect on plasma jaotion in the stationary heating regimes 
with one exception, the case of a picosecond pulse with a CO^ 
gas laser. 

a. Calculation of Electromagnetic Forces 
The dispersion relation for electromagnetic radiation 
subject to inverse-hremsstrahlung abosrption is. 


n = n - ih = [l - 


-f- 1/ T 


1 

HIT 


, 1/2 


+ 1/ T 


where n is the complex index of refraction, tu and are 

the laser and plasma frequencies respectively, and T is the 
electron-ion collision time. For a sufficiently slow variation 
in the index of refraction, the approximation of geometrical 
optics can be used.^ 

The Loren tz force density can be written in terms of the 
Maxwell stress tensor and the Poynting vector. Following Hora, 
this force is. 


? * 1 , „2 2 _ 1 _ 
= - x 1^2 (5E +uH > “ c 2 


bt 


E x H 


for a one-dimensional plasma with a density variation in the 
direction with a plane electromagnetic wave traveling in the 


x 
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+ x direction. 

Only the time averaged value {over a few oscillations) 
of the force density is considered. Also, the time average 
of the Poynting term will be negligible for laser pulses 
much longer than the period of the wave. Then the time 
averaged force is. 


f B i ISl-S- e -*r 

em X 1 4 ax e 


{Re n“ 1/2 ) 2 
(Re n 1/2 ) 2 + (im n _1/2 ) 2 + (1m n 1/2 ) 2 ] , 


where k represents the damping of the wave, 

/ 

/ , x 

k = f K dx 
j c J 0 

i 

and .Ea is the amplitude of the electric field of the radf 
ation entering the plasma. 1 /ujt will be small for plasmas 
irradiated with a laser beam in the visible or longer wave 
lengths for a temperature greater than about 20 eV. The 
amplitude of the electric field is simply related to the ' 
intensity — and to the terminology introduced in section A, 




r 


where = d$/dT . Moreover, the solution to the radiative 
transfer equation {10} is, 

_ i « S'e”* 
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Hence, the magnitude of the time averaged electromagnetic 


force is. 


em 


J. /A , r 

0 f|i(l/n + n) + 


2Rct ayL 



where 

n = (1-C) 1/2 + 0[1 /(ujt) 2 ] (2-77) 


In the stationary- frozen-conducting regime, in which 

i = § 1 [1 + 0 (a 2/5 ) 3 , 


f..„ = 
em 


J 0 A 

4Rct 


( 1-0 




3/2 dy 


(2-78 ) 


The apparent singularity at the critical density is not a real 
one since the expression used for the index of refraction (2-77 } 

fails near that point. The stationary-frozen-conducting regime 

is presented since it is in that regime, if at all, that elec- 
tromagnetic forces may he important. 

b. Calculation of Thermal Forces 
The thermal force density is simply the pressure gradient 
Vp . Assuming the contribution of the cold ions to the pres- 
sure is negligible, the thermal force is. 


M 


<n e *T e ) = 


J o /A 

R 2 


_b_ 

Oy 


(C8 e ) 


In the stationary- frozen-conducting solution from (2-52) , 



«!«> 


2/5 &L 

dy 
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■where the ratio of the two integrals in (2-52) was left out 
for brevity. This ratio is of order one if the dark part of 
the plasma is roughly the same size’ as the light part, 
c. Comparison of Thermal and Electromagnetic Forces 
The ratio of the forces is. 


- [ E/c yf a)2/5 ] 7575 ltvs 


•-ft 


C* - ' ''(1-0 3/10 (7-2£) J 


(2-79) 


2/5 

The first part of the expression, R/ct {5/2 a) ' , gives the 

2/5 

order of magnitude of the ratio. The second part, $'/$ 
gives the time, variation and is generally of order one or less. 
The third part, a complex function of £ , gives the spatial 
variation, and it is of order one except near the critical 
density, £ — 1 , where there is an apparent singularity. This 
is not a real singularity of course since the expressions used 
to derive (2-79) breakdown near the critical density. It is 


clearly seen that — holding other things constant — it is possi- 
ble to make ~f /Ivpl as large as one wants if the pulse time 
t is made small enough. 

The first factor is in the usual terminology. 


R/ct (fa) 2/5 =1.15 X lO'^J/A) 375 !^" 1 


Given X . t , then when J„A is larger than a certain 

0 * p ' v 

value, the factor will exceed one- This gives the approximate 

size of J A for which electromagnetic forces exceed thermal 
0 

forces. The table below shows these values of J Q A 
typical pulses and wave lengths. 
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1 nsec 

1 psec 

i.06 p 

13 

1.6 x 10 

1.6 x 10 8 

10.6 p 

9 • i 
7.4 x 10* 

7.4 x 10 4 
2 



joules/m 


Clearly then, the electromagnetic forces will exceed the thermal 

forces for values of J A well within the range of interest 

0 

to laser plasma heating. 

It is established then that electromagnetic forces may 
well exceed thermal forces. But this does not mean that they 
will have a major effect on the plasma motion, especially 
since the largest electromagnetic forces occur in the pulses 
of shortest duration. The real measure of their effect must 
be based on the total motion they generate. 

d. The Overall Effect of Electromagnetic Forces 
y;hen a force acts on an object for a certain length of - 
time, the object attains a certain velocity. Then the effect 
of electromagnetic forces is to impart a certain velocity to 
individual particles of the plasma. If that velocity rs small 
compared to the local acoustic speed, a , then the effect of 
electromagnetic forces on the subsequent motion is small. If 
that velocity is of the’ same order or larger than a , then 
the electromagnetic forces strongly affect the motion. 



. m 

. A calculation is made of the velocity, V em , generated 
in a nonuniform plasma by electromagnetic forces. For simpli- 
city, it is assumed that only electromagnetic forces operate, 
and that each individual particle is accelerated by the force 

for the time t . The effect of any motion during the laser 
P 

pulse on f is neglected. The calculation uses as the 
equation of motion, 

dV 

P = f (2-80) 

K dt em 

Integrating (2-80) from t = 0 to t = t using f em from 
(2-78) and an acoustic speed based on the stationary-frozen- 
conducting solution, the ratio of the electromagnetic to the 
acoustic velocity is found. 


4ra.cn Ra (| a) 2 5 
i ec o ^ 


(1-C) 7/ V /5 


(2-Ql) 


where a Q - is the acoustic speed based on the temperature T e . 

The size of V /a is roughly the size of the term in brackets 
em 

and is given in figures 31 and 32 for a Nd glass laser and a 
C0 2 gas laser respectively. 

It is recognized that the spatial variation given by the 
term in parenthesis may be much larger than one in two places, 
near the plasma edge, £ = 0 , and near the critical density 
C = 1 . The motion near £ = 0 will be dominated by the dyna- 
mics of the adjacent denser parts of the plasma. The electro- 
magnetic force does seem to have a significant effect near 
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q - i though, but the presence of anomalous absorption effects 
in this neighborhood will likewise generate high thermal forces 
which may well be more important. The anomalous absorption 
effect has been studied theoretically by Kriier et al.' L£ ’ and 
was seen to appear in an earlier experiment by Gekker and 
Sizukhin.^" 7 

The comparison of figures 31 and 32 with figures 16, 17, 

18, and 19, the re-examined regime charts shows the regimes 

in which electromagnetic forces become important. It is seen 

that V /a becomes of order one only in a heating regime 
em 

that is already nonstationary and hence not of interest to 
the problem of stationary laser plasma heating. The one excep- 
tion is a picosecond pulse with a COj laser where electromag- 
netic forces become important at the edge of the stationary- 
frozen-conducting regime. But even then the nonstationary 
regime is only extended slightly. 

Appendix 2 

The Lawson Criterion for Heating a Uniform Plasma 

The focus of laser plasma heating as described in this 
work is to achieve a thermonuclear reaction that releases more 
energy than required to generate it, i.e., the achievement of 
gain. Thus it is appropriate to conclude this chapter with 
pertinent statements regarding the achievement of gain. 

The standard requirement for gain in a pulsed reactor is 
known as the Lawson criterion.^ The. Lawson criterion is that 
the product of the number density of the reactants (equal 
parts Deuterium and Tritium) and the confinement time must 
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exceed 10 14 sec/cm 3 , and that the temperature of the reactants 
exceed 10 KeV, 

nt > 10 14 sec/cm 3 (2-82) 

T. > 10 KeV (2-83) 

x 

Physically this means that a sufficiently dense D-T plasma 
must react for a long enough time at a high enough temperature 
in order to release more energy than went into heating the 
plasma. The beauty of the Lawson criterion is seen by its very 
simplicity — and by its applicability to most any scheme or 
geometry. Of course, when the design stage is approached, a 
more exacting criterion will be needed to take into account 
the details of the thermonuclear reactor. For the analysis of 
a laser heated plasma the confinement time t is taken to 
be the ordinary acoustic time, , from (2-76) . The number 

density of the reactants is " n e ~ x 10 //cia 

Thus the first part, (2-82) of Lawson's criterion becomes, 

j /h > 3.60 x 10 12 ( ; — 
o' v 

The second part, (2-83), of Lawson's criterion is that the ion 
temperature exceed 10 KeV. This is true if the electron tem- 
perature exceeds 10 KeV, and if there is electron-ion equili- 
brium t > t , (2-70)'. These are presented in figure 33 
where the various requirements are. plotted against laser 

j /a , and the density, Q . 


energy per unit area. 



Joules/m' 


V A 


Figure 33. Lawson's criterion for a 
uniform laser heated Dlasma 
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The first part of Lawson's criterion requires £ to be 
14 

below the line nt = 10 . The second part requires q to 

be below the line t = t and above the line T > 10 KeV. 

Lines are also shown supposing that the confinement time is 
ten acoustic times. 

11 2 

The minimum energy for achieving gain is 10 joules/m , 
large but probably feasible for future lasers. The corres- 
ponding density is 0.1n ec and the plasma length is 60 m. For 

18 —3 

a CO,. gas laser this corresponds to n = 10 cm .A mag- 

netic field of about a megagauss would be required to contain 

such a plasma. This is a large magnetic field, but not 
beyond possibility for the future. A 60 m length for such 
a magnetic field is a major engineering project to be sure. 

If it were possible by some technique to confine the 
plasma for ten acoustic times (by a magnetic mirror effect 
possibly) , then the minimum energy needed reduces to about 
10 10 Joules/m 2 , corresponding to a density of 5 x 10 18 for a 
CO^ gas laser. The confining magnetic field needed is only 
increased by a factor of two, and the plasma length is a very 
reasonable 60- cm. This is an attractive possibility if the 
2 MG magnetic field together with the enhanced confinement 
can be -attained. 


Appendix 3 

Effect of Z-Doping on Lawson's Criterion 
One of the characteristics of laser heating with a long 
wave length laser (such as the 10. 6p CO 2 gas laser) is a 
fairly long absorption length. A long absorption length 
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requires that the plasma volume be large and thus the laser 
energy required will also be large. 

A possible solution to the problem of unwieldy absorption 
lengths has been suggested by Dawson, 22 and that is to dope 
the D-T plasma with an impurity having a high atomic number. 

The presence of the “high Z" impurity will increase the elec- 
trons significantly and thus increase the electron-ion collision 
time t . Increasing t shortens the absorption length. 

2/5 

Specifically, the absorption length is proportional to 1/Z 
where Z is the effective atomic number of the ions. Z = 1 
for a pure D-T plasma. It is appropriate then to examine the 
Z-doping approach in the light of achieving gain in a thermo- 
nuclear reaction. 

First of all, the temperature will go up — simply because 
a shorter length of plasma is being heated. Raising the tem- 
perature will also raise the pressure. Now if the plasma is 
inertially confined, then an increase of pressure may be of 
no consequence.. But in a magnetically confined plasma, only 
so much pressure can be contained. Thus the plasma must be 
initially at a lower density so that the higher pressure can 
be contained by the magnetic field. But reducing the density 
to bleep the pressure fixed increases the absorption length 
and the net effect is that the absorption length is not 
reduced at all! Thus, Z-doping is of no value in a magnetically 

confined plasma. 

Subsequent discussion assumes that the electron number 
density is held constant since it determines the refractive - 
properties of the plasma (when reflection occurs, etc.) . Then 



119 


the length of the plasma is reduced. This will immediately 
reduce the confinement time, which is detrimental with respect 
to achieving gain. Moreover, doping with an impurity (holding 
constant) reduces the number density of the reactants — 
and hence the thermonuclear energy release. Also, bremsstrah— 
lung losses will increase due to the increased ratio of elec- 
trons to ions. 

A measure of the effect of Z-doping on gain has been made. 
There is a certain fixed amount of energy that must be added 
to a pure D— T plasma of a given electron density to raise its 
temperature to 10 KeV. A minimum nt is then needed to 
achieve gain from the reaction that follows. Suppose that 
exactly .the same amount of energy is added to a doped D-T 
plasma which has the same electron density. Then a 
new minimum nt for the doped plasma can be calculated using 
the D-T reaction and the bremsstrahlung power densities. 



br 


= 2.59' x 10" 24 ^§73 f exp - ^t7T ] a l wafct / cm3 
= 5.35 x 10“ 31 n 2 T ly/2 . a 2 watt/cm 3 


where n is the number density in /cm 3 of all particles, and 
T is the temperature in KeV. An equilibrium plasma is assumed, 
a. and a- are factors depending on the doping, 

i. * 

ci ± f(R-i)in + f(A 2 -m 

3l [1 + 5*f (A-l) ] 2 
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L = . U 

2 [1 + ?2f(A-l)] 2 


f is the impurity fraction of the ions and A is the atomic 
number of the impurity {assumed fully ionized). Then Lawson's 
criterion is 

_ J-|F- (2-84 ) 

nt 1-tj 3n \ T 3n 2 kT 


r| is the efficiency of conversion of released energy into 
heating the next slug of plasma to thermonuclear temperatures . 
This calculation takes g = 1/3 after Lawson. Using the 
above equations, the required n g t for gain (i.e., satisfying 
the inequality . (2-84) was calculated and is plotted m figure 
34 . Used in the calculation was n = a^n e where 

" ‘ 1 ~ f 

a 3 " 1 + f(A-l) 


Also used was T = a.T where T is the temperature that is 

*+ o 

achieved in the doped plasma with the same amount of energy 
that produces the temperature in the pure plasma. a 4 

is given by 


4 


Cl ± f (A-l) 3 7/ 5 

tl + %£ (A-l) J 


It is clearly seen that Z-doping always makes the Lawson 
criterion more difficult to attain , i.e., the minimum n g t 
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ion fraction of Impurity f 

Figure 34. Lawson criterion for a high-Z doped DT plasma 
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aj^ays raised. Moreover, Z— doping does not improve the 
electron-ion equilibration time. In fact it raises t g( ^ 
(2-76} , since is increased, and n g is held constant. 

Thus the conclusion is obvious; Z-doping is to be discarded 
insofar as the goal of a thermonuclear reactor is concerned 



Chapter Three 

LASER REFRACTION IN CYLINDRICAL PLASMAS 
A, REFRACTION IN A STATIONARY PLASMA COLUMN 
I. Introduction 

One of the' promising geometries for the production of a 
thermonuclear reaction is a long cylindrical plasma column. 

Such a column might be contained on the sides by a magnetic 
field which could be generated in several ways. But magnetic 
fields can only hold so much pressure, thus the plasma density- 
must be made low enough so that the pressure at 10 KeV can be 
contained. Previously this -was an impossible drawback to the 
laser heating of magnetically confined plasmas since the cur- 
rently available magnetic fields (ca. 300 KG) could only con- 

17 3 

tain a plasma density of ca* 10 /cm * The absorption length 

17 3 

- for a Nd+ glass laser in a 10 /cm , 10 KeV plasma is an 
incredible 450 miles l 

However, the advent of powerful long wave length lasers 

has made a dramatic change in this situation. Dawson, Kidder, 
23 

and Hertzberg have proposed the use of a long wave length 

gas laser to heat such a magnetically confined plasma. Recent 

24 25 26 

developments of high power pulsed CO^ lasers 1 * with 

.lo = 10. 6u show that gas laser development is becoming compet- 
itive with the solid state lasers in the production of high 
energy pulses. The significant advantage in using the long 
wave length, as first noted by Dawson, is that the laser energy 
can bo better absorbed by plasmas at densities which can he 
confined in the magnetic fields that are currently feasible. 
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The mechanisms of absorption (with the exception of the 

16 

recently discovered anomalous absorption 3 are fairly well 

understood. The primary mechanism is inverse bremsstrah- 

l-ung^ 4 '^ 7 ’ In addition, nonlinear^® and resonance phenora- 
11 29 

ena ’ have been examined. The actual character of "the 
heating has also been studied (see references cited in Chapter 
II). But with one exception, nothing has been said about 
refractive effects in a columnar plasma: Vlases and Ahlstrom 30 

studied the thermodynamic, absorption, and -refraction problems 
associated with heating a magnetically confined plasma column 
with a long wave length laser. 

’ This lack of previous work is not to indicate that refrac- 
tion is not a problem. On the contrary, an attempt has been 
made to propagate a laser down a typical magnetically confined 
plasma (ft -pinch) . The result was that the beam was refracted 
out of the column to the extent that no signal penetrated to 
the far end of the plasma . Thus the refraction may play a 
very significant role in plasmas of this geometry. If it is 
not possible to keep the laser beam in the plasma column, then 
all the advantages of long wave length heating are lost, and 
it is critically important to study the refraction problem. 

' - a. General Refraction Problem 

When the geometry of a particular scheme is investigated 
it becomes obvious that refraction effects can play a dominant 
role. The density gradients in the plasma may refract the 
radiation out of the plasma before it has traveled far enough 
to be absorbed. The equation for the index of refraction in a 
gas composed of atoms, ions, and electrons is 
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n = 1 + (Po+Ko '\o 2 )n c + l K i +K i , *o 2 ) n i - 4.5 x 10~ 22 n e X 0 2 , 

-3 

where n 0 « cm and Xo ~ microns. This equation 
applies when the frequency u> is much greater than the plasma 
frequency uj^ and the cyclotron frequency. For x© > 0.5y 
and a typical gas the contributions due to the atoms and the 
ions can be combined so that 

n « 1 + lO^tno+rVj.) - 4.5 x 10~ 22 n e x o 2 . 

Thus it is clear that for Xo > 0-5u and ionization greater 
than 10?S that n - 1 < 0 whereas for an un-ionized gas 
n - 1 > 0 . Therefore in an un-ionized gas the refraction 
deflects the beam towards the higher density but in a fully 
ionized plasma the index of refraction decreases with increas- 
ing density so the refraction deflects the beam towards the 
lower density. Since all plasmas have n g -» 0 at their 
boundaries, refraction could play a dominant role in the 
application of laser radiation to the production of CTR. 

The simple equation relating the radius of curvature of 
a ray, R , to the index of refraction gradient is 

1 _ sin ft i . ~ 1 
n ~ ~ 1 9 ra< ^ n 1 t 

** n 

where ft is the angle between the ray and the direction of 
grad n , It is particularly interesting to note’ that r o 
at the plasma critical density, i.e., where the laser fre- 
quency equals the plasma frequency. For laser heated plasmas, 
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the maximum density P roax , may be greater or equal -to the 
critical density, p c . 

Jn this study only the problem of a cylindrical plasma 
with azimuthal symmetry irradiated in the axial direction is 
discussed. The results developed here indicate effects that 
would be important in other geometries as well, and the general 
trends can. be carried over. 

b. Unfavorable Density Profile 
The most obvious geometry, where p{r) is parabolic with 
Pmax at b = 0 , is considered first. The refraction effects 
are detrimental because of the unfavorable density gradient. 

It is shown that heating this type of plasma is very difficult 
due to the refraction losses. 

If the incident laser beam is focused then it is shewn 
that there exists a solution Where the laser beam can be made 
to propagate indefinitely even in an unfavorable density 
gradient.’ Realistically some deviation from the ideal case 
would occur in practice. When a realistic deviation is allowed 
then it is shown that the propagation distance is improved by 
a factor of approximately three over the unfocused case. 

c. Favorable Density Profile 
The other possibility which exists is that the laser beam 
is propagated in a region where the density gradient is favor- 
able, i.e., leads to self focusing. It is shown that for a 
favorable density variation a parallel incident beam is trapped 
as in a light pipe. Solutions for the wave fronts in a para- _ 
bolic density profile are presented. These solutions based on 
geometrical optics show that the wave fronts are a family of 
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cusped curves that extend down. the column, alternately collaps- 
ing and expanding. During this process envelopes of the wave 
fronts, caustics, are formed at the locus of points of the end 
of the cusps. 

It is also shown that m a light pipe configuration of the 
plasma, the effective absorption length of the plasma can be 
reduced by as much as an order of magnitude. 


2. Defocusing in a Parabolic Density Profile 
a. Unfocused Laser Beam 

The geometry considered is a long cylindrical plasma 
column with azimuthal symmetry and a density maximum on the 
axis with a monotonic decrease of the density in the radial 
direction. A laser beam is incident on the plasma column in 
the axial direction where it is assumed that 'the diffraction 
effects can be neglected and the laser beam is in. the TEM^ 
mode with the intensity distribution approximated. by a para- 
bola. Phe general configuration is shown in figure 35. It 
is accurate to use geometrical optics as long as the plasma 
density is less than the critical density and caustics are 
not generated. The equations for the intensity and the 
abosrption must be modified near the critical density and 
caustics if the details are desired in these localized regions. 

The Eihonal equation for a cylindrically symmetric geome- 
try is 


c 2 {r) 


( 3u 
V bx 


) 2 



where u corresponds to time and c is the phase velocity. 
This equation can be reduced to a differential equation for the 
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dx 

where the two boundary conditions are that 

r(0) = r 0 , 


;( 0 ) = 0 . 


The integration of (3-1) leads to the ray equation 


f r car 


The phase velocity is given by 


c « Cq (1-C) ^ , 


where c, is the speed of light in a vacuum and < = ^2 " P c ' 
the ratio of the density of the plasma to the critical density 
for the particular laser wave length. Thus (3-3) becomes 

r\ — T — r r . (3-4) 

* * ^ J r 0 • 

where Co = C ( r o ) and C = C (*> ■ 

Now in order to evaluate the integral in (3-4) , a para- 
bolic density variation is assumed. 


c - cl 1 - r i> 2 l • 


(3-5) 


T * i corresponds to C = 0 at the edge of the plasma r - b; 
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and r < 1 corresponds to C = (1-rj C ra at • r = b . Thus for 
r < 1 the density profile is flatter as shown in figure 36. 
The ray equation, (3-4), using (3-5) reduces to . 


r 


cosh 


M— 


rr 


(1 -T r 0 /b) 


■i) 


(3-6) 


(3-6) can then be used to determine the distance that the 
laser beam propagates along the column. A measure of this 
distance is the value of x for which one half the laser beam 
energy has left the plasma column, x^ . Figure 37 shows the 
variation of Xj^/b with T for three values of C m * F° r 
each of the values of C m a value of the absorption length 
for X 0 = 10. 6 u and = 10 KeV is also given. It is clear 
that the radius of the plasma .would have to be very large for 
any significant portion of the energy to he absorbed at 10 KeV 
and that F should be very small. 

It should also be pointed out that the case F < 1 could 
be considered as a full parabolic plasma profile, T = 1 , with 
a smaller radius laser beam r fa . For example if b is the 
plasma radius then r^/b = "3 corresponds to T - Yj-j and 
r b /b = Yfi corresponds to r = Yqq • As lon 3 as the character- 
istic thermal conduction time in the lateral direction is much 
less than the longitudinal acoustic time, the whole plasma 
column would be heated. 

Finally for this case it is straightforward to roafce a per- 
turbation calculation of the percentage of the energy absorbed 
at any given temperature. Figure 38 shows that for a CO ^ laser 
and temperatures greater than 100 eV a significant fraction of 
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Figure 37. Distance in which half of the laser energy leaves 
the plasma, for a parabolic density, and intensity profile 


13 3 



Figure 38. Energy absorbed in a parabolic density 
■ profile at a given temperature 
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the energy is not absorbed and at 10 KeV the percentage 
absorbed is negligible. These calculations were done for 
b = 1.0 cm and T = 0.10 . For larger plasma diameters and 
smaller values of V the picture would of course be more 
favorable. 

b. Focused Laser Beam 

Now suppose for the parabolic density profile, the 
incident laser beam is focused by a lens such that in the 
absence of the plasma the focal point would be at x = h as 
shown in figure 39 - Before solving this case two limits are 
obvious. If | - - then this is 'the unfocused case dis- 

. . t* b n then the radiation 

cussed in the previous section. If b 

will traverse the plasma column and go out the other side. 
Solutions for finite values of h/b should give increased 
propagation distances , in fact it would seem reasonable to 
search for a solution such that the radiation is trapped. 

The boundary conditions are n.ow 


r (0) - r o t 


dr 

dx 


to) = - 


(3-7) 


The 


solution for this case has three forms: 


{ - vtsinlT 1 sfr - sintT 1 bftf £ ° r 9 " ° ' 


$ = Ytlnf - in §•} 


for p - 0 , 


(3-8) 

(3-9) 




f = Ytcosh" 1 - cosh' 1 sflj 1 


tor p < 0 , 


(3-10) 


where 


p 2 =. 


.1 + 


h‘ 


"7^ 

m b 


and 




1 -5- 


_ 2 1 r ' c« 

To m 


1 ) + ^1 
b 


For p 2 > 0 the ray paths are overfocused so that they cross 
over the axis and exit from the plasma on the other side. For 
g 2 < 0 the ray paths are underfocused (as shown in figure 35) . 
For p 2 = 0 all the ray paths approach the axis asymptoti- 
cally. This condition requires 


h_ = i ( i _ 1 ) 

,2 r 1 C 1 


(3-11) 


ra 


It should be noted that in heating a plasma, this type of solu- 
tion would only he useful if the characteristic thermal conduc- 
tion time in the lateral direction is much less than the longi- 
tudinal acoustic time so that the whole plasma is heated. 

As a practical matter it would be impossible to satisfy 
(3-11) exactly. Thus it is of interest to examine the effect 


of an error in e.ct. 


and (3-10) give 
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where r c is the correct value for (3-11) to be satisfied. 

It is easily seen that only a small error in specifying r 

Xi 

leads to values close to those obtained in the unfocused 

case. Figure 40 shows the variation of vs T for r - T c = 

0 .lor c and 0.0ir c . It is seen that Xj^/b is improved 

over the unfocused case by roughly factors of three and nine 

for I* - r = o.ior and O.Oir respectively, 
c c c 

3. Refraction in a Favorable Density Profile 
It is clear that if the density profile has a minimum the 
rays will be refracted into this region. It may even be pos- 
sible to trap the rays in a minimum density region. The con- 
figuration considered is the same as that of Section A 2 except 
now it is assumed that a density minimum exists. In this case 
it is both convenient and interesting to use the equations in 
their characteristic form so that 


and 


dr 

ds 


+ 


2 


C( r 0 ) - C(r) 



> 


dx _ . 2 
ds c 0 


( 1 - C(r 0 ) ) . 


(3-12) 


(3-13) 


Consider the density variation shown in figure 41 . From 
(3-12) it is clear that any ray entering the column with 
r s: r 0 s r 2 and ||(0) = 0 will be trapped in this minimum 
region. Such a configuration could then be called a "light 
pipe" since the light is trapped in the plasma column. 

In order to study the behavior of the radiation in a 
minimum region in more detail, it is instructive to assume a 
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Fiirure in. Density Variation Kith a general mininma 











s and t are two parametric variables. The wave fronts arise 
for constant values of the time variable, u . 

The general proof has already shown that all the rays v.culd 
be trapped and is confirmed by (3-15) . For this parabolic 
density profile all the rays are cosine functions where the wave 
length of the oscillating rays goes to zero as r 0 -* b . 
Physically, this means that those rays entering the column near 
the critical density point are bent the most sharply and 
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oscillate back and forth across- the column with a very short 
wave length. Those rays that enter the column near the center 
are only weakly deflected and oscillate hack and forth very 
slowly. The ray that enters the column right on the axis will 
proceed along the axis without oscillation. 

A series of wave fronts are shown in figure 42 . They are 
a series of cuspad curves that alternately collapse to the 
axis and then expand. This behavior leads to a series of 
envelopes or caustics at the points of reflection, i.e., where 
the wave front turns around and heads back toward the axis. 
Several caustics are shown in figure 43 . 

These solutions have infinities on the caustics which are 
due to the geometrical optics approximation. The behavior 
near a caustic is discussed in the first appendix. The axial 
symmetry also leads to a singularity at the axis. The behavior 
near the axis is also discussed in appendix one. 

4. Absorption Length with a Favorable Density Profile 
Xt is of interest to compare the. absorption length in a * 
uniform plasma to the absorption length in' the plasma with a 
favorable density profile. Naturally, the absorption length 
in the "light pipe" is not a clearly defined quantity. The 
rays of the laser beam that enter the central region where the 
density is lowest will experience weaker absorption than those 
rays which enter near the periphery— at a higher density. 
Furthermore, all rays which do not enter at the density minimum 
will oscillate back and forth across the plasma as they move 
down its length. Thus, at any point in the column one finds a 
strange conglomeration of rays, whose intensity has been damped 
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Figure kZ. Wave fronts in a plasma column 
with a favorable parabolic density profile 
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t>y widely varying amounts, and whose position seems to have no 
relationship to the origin of the rays. 

The appropriate way to measure an "effective absorption 
.length for such a plasma is to calculate how the overall laser 
power diminishes in moving down the plasma column, lie approach 
is to find first how each bundle of rays has diminished in 
intensity. Consider an annular bundle of rays of equal inten- 
sity, 

dW(x,r 0 ) = l(r 0 ) * 2Trr 0 dr 0 exp|-j^K[r(x,r 0 ) ]ds} _ 


Where dw(x,r 0 ) is the power at x of the radiation that 
entered the plasma between the circles r = r 0 and 
r = r 0 + dr 0 . I(r 0 ) is the intensity profile of the laser? 

K(r) is' the absorption coefficient, and ds is an incremental 
distance along the ray which began at r 0 . The total power 
in the cross section at x is found by integrating dw over 

all r 0 > 

■ - • b r x i 

w(x) = 2 tt J 1 I (r c ) * exp^-j^ Kdsjr 0 dr 0 . (3-17) 

The effective absorption length can be calculated by performing 
the 'integral (3-17) for a particular density profile. Weak 
oscillatory terms will arise in the calculation of W(x) but 
these should be neglected giving W(x) . The effective absorp- 
tion length is then given by 


'- off = (- w t x >/^^)x=C 


(3-18) 
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Calculation of ^ e f f^ 0 T ei was P er ^° :nnei ^ f or t ^ le favorable 
parabolic density profile (3-14) using two different intensity 
profiles? I{r) ~ 1 - (r/b) 2 which approximates the TEM Q0 
mode, and I{r) ~ (r/b) [1- (r/b) ] which approximates the 
TEMqq + q^ mode. c 0 is the speed of light in a vacuum and 
is the electron-ion collision time. The results are presented 
in figure 44 where they are compared to the absorption length 
in a uniform plasma at a density equal to the minimum density 
of the favorable profile. It is seen that there is significant 
enhancement of the absorption length. In fact there is an 
actual maximum absorption length which occurs approximately 
when the minimum density is half the critical density. There 
i 3 at least an order of magnitude decrease in the absorption 
length for all densities less than about 0.4 of the critical 
density. 

Thus, it is seen that the light pipe not only traps the 
laser radiation, but reduces the absorption lengths considerably 
from the often inconveniently long lengths which arise in a 
uniform density plasma. 


5 . Summary 

The problem of propagating a laser beam along a long cylin- 
drically symmetrical column of fully ionized plasmas has been 
considered using geometrical optics. For the case where the 
density has a maximum on the axis and decreases monotonically 
in the radial direction the ray paths are refracted out of the 
plasma column. The solutions show that it is very difficult 
to heat a long plasma column with this unfavorable density 
profile. By focusing the incident laser beam it is theoretically 
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possible to keep the light from being refracted out of .the 
plasma column. However, this solution can be compared to 
finding a neutral stability point in the middle of an unstable 
region. So that any deviation leads to a very rapid refraction 
Of the laser beam out of the plasma . 

A "light pipe" effect is found if there is a density- 
minimum. For the general case of a density minimum it is 
shown that a parallel incident beam is trapped by the plasma 
column. A special case where the density variation is para- 
bolic is considered, and the solution shows that the wave 
fronts are a family of cuspea curves which alternately collapse 
to the axis and then expand until reflected back along a caustic. 
The intensity along the caustics is very large as is the inten- 
sity along the axis. Geometrical optics fails in these local- 
ized regions (see appendix one) . 

Finally it is shown that for the “light pipe" solutions 
there is a significant decrease in the effective absorption 
length due to the rays oscillating back and forth across the 
plasma. It is clear from these solutions that there are very 
significant advantages to the' light pipe configuration for the 
production of a thermonuclear plasma. 
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B-- DYNAMIC GENERATION OF FAVORABLE DENSITY PROFILES 
1 . Introduction 

Since the plasma density profiles generally encountered 
are unfavorable with respect to refraction (such as in the 
quasi-steady stage of a 0-pinch) , it is important to look for 
ways of creating a favorable profile. There are basically 
two approaches to this problem. ‘The first approach is to catch 
the plasma during its formative stage when the density may not 
have settled dovm to its steady "unfavorable" state. The second 
approach is to take the quasi-steady unfavorable profile and 
change it to a favorable one, i.e., density profile tailoring. 

The first approach can be applied to a pinch device. Dur- 
ing the formation of the pinch, an axially symmetric shock wave 
collapses to 'the axis. Across the shock the density rises. So 
at any instant of time during the collapse, a favorable density 
profile is seen to appears a lower density on the axis than 
farther out. The density profile during the "collapse phase 
would resemble the solid line in figure 45. 

After the collapse phase comes the "bounce phase" when the 
shock wave reflects from the axis. ‘Due to the nature of the gas 
dynamics, a favorable density profile is left behind the reflec- 
ted shock. The density profile corresponding to the bounce 
phase is the dashed line in figure 45. 

If laser heating during the collapse phase is used, the 
laser pulse must be shorter than the time for the imploding 
shock to collapse through a distance equal to the radius of 
the beam. Likewise, for the bounce phase, the pulse must be 
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Figure ^5. Favorable density profiles In the collapse 
and reflection phases of a 0-pinch 
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shorter than the lifetime of the favorable density profile. 

All subsequent discussion will deal with the second approach 
that of changing an unfavorable density profile into a favor- 
able one. The basic method here might be called "selective 
heating" whereby the region near the axis is heated with a small 
laser, the “pre-laser" . The added energy would produce a, pres- 
sure imbalance causing redistribution of the plasma and possi- 
bly a favorable density profile for trapping the beam from the 
main laser. 

This approach offers a great deal of flexibility in achiev- 
ing the desired' result. The pulse lengths, the wave lengths, 
and the energy of each pulse can be treated as independent 
parameters. For example, suppose the intent is to prepare a 
plasma for heating by a powerful CO 2 laser. Then the plasma 
length would correspond roughly to the absorption length of 
10. 6p radiation at 10 KeV. But the absorption length at, 
say^ 100 eV is only 1/1000 of its value at 10 KeV. Hence, a 
CO 2 pre-laser could not be used to heat the plasma longitu- 
dinally- Only a small fraction of the length would get heated. 
But a Nd + glass laser which has a longer absorption length may 
be used. 

Another advantage of using a shorter wavelength for the 

pre-laser is the refractive improvement. Figure 37 shows that 

if the initial plasma column has n e x = 5 x 10 /cm and 

r ** 0.10 , then for a C0 2 gas laser, one half the laser beam 

would have left the column in three plasma radii.. However if 

* ^ 

X = 1.06p corresponding to a Nd + glass laser, then x^/b ~ 100 
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Figure 46 (which is the same as figure 38 except X Q = l.OSp , 

Nd + glass laser) shows that if b - 1 cm , 8056 of the energy 
would be absorbed. 

a second consideration is that the pre-laser preparation 
of a favorable density profile need not be a high efficiency 
process. For example, if the initial temperature is 100 eV , 
then this represents only 1% of the energy required to achieve 
10 KeV. If the heating from 100 eV to 200 eV were done in a 
way that produces a favorable density profile, then the remain- 
der of the heating (to 10 KeV) could be done very efficiently, 
since there would be no refraction losses. Inefficiencies in 
the heating to 200 eV could be tolerated since this represents 

only 1% of the total energy required. 

It is also interesting to note that once a favorable den- 
sity profile has been produced the additional energy absorp- 
tion due to an incident Gaussian las'- r beam will continue to 
push^mass towards the outer edge of ie beam. This problem 
is very similar to the thermal blooming problem 31 ' 3 which is 
experienced in the propagation of a laser beam through the 
atmosphere. The fundamental and convenient difference is that 
in the case of a fully ionized plasma the laser beam is trapped 

because n - 1 ~ -n g • 

2. Governing Parameters 

Before considering the methods of density tailoring, it 
is necessary to calculate the various parameters that govern 
the dynamics. Rather than set up the full equation, it is 
sufficient at this point to examine the scale- of the parameters. 
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The length scale is 2r fa where r b is the radius of the laser 
beam entering the plasma, r b will be less than or equal to 
the plasma radius, b . The density and temperature scales 
are the maximum density and temperature respectively in the 
unheated plasma column. 

There are several time scales of interest. There is the 
heating time scale t^ which is the time in which the tempera- 
ture would double if the plasma were held stationary. Then 

T „ , 1 K v W 
** = C v P nr b 2 ' 

where T is the initial temperature, C v is the constant 
volume specific heat, K v is the absorption coefficient, p 
is the density, and W is the laser power. Then t h becomes 

^ = 7.40 x 10 17 r b 2 T 5/2 / W\ 2 n e sec . {3-19) 

where r b is in cm , T in eV , X g in microns , W in 
, -3 

watts , and n e in cm , 

The acoustic time scale t Q is defined in the usual way 
as the quotient of the length scale and the acoustxc speed. 

Then 

fc a = 1.80 x 10“ 6 r b t“ X 7 2 sec (3-20) 

The thermal conduction time t^ is defined as the time 
for significant thermal conduction to traverse a distance 2r b , 

fc th “ pC v t2r b )2 7 K e > 
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where K is the electron thermal conduction coefficient, 
e 

Then 

t.. = 1.29 x 10~ 20 r. 2 n t' 5/2 sec (3-21) 

t:n o e 

Figures 47 and 48 compare the acoustic and thermal conduc- 
tion time scales to the heating time for \ Q - 1.06p and 
10. 6n , respectively. Figure 49 compares the acoustic time to 
the thermal conduction time. The interpretation of the various 
limits will depend on the particular method of density tailor- 
ing, but a general comment can be made. Xf the thermal con- 
duction time is shorter than the characteristic time of the 
particular method, then thermal conduction will tend to wipe 
out any favorable density profile that might be generated. 

This conducting regime poses a major difficulty. A nonconduc- 
ting regime occurs if the thermal conduction time is longer 
than the characteristic time of the particular method, then 
favorable density profiles generated will endure at least for 
that characteristic time. Several schemes for density tailor- 
ing are presented below for both conducting and nonconducting 
regimes. 

One more time scale must be introduced, the magnetic 
diffusion time, which will be important in determining the 
characteristics of the conducting - ^gime. The magnetic diffu- 
sion time, tg , is the time in which a significant amount 
of the magnetic field diffuses through one characteristic 
length. It is found by an examination of Maxwell's equations, 

" ^< 2 r b > 2 
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Figure 47. Comparison of acoustic and thermal conduction 
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time scales with the heating time scale for a Nd glass laser 
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where p is the permeability, and o is the electrical con- 

14 

ductivity. Using electrical conductivity from Spit2er, fcg 
becomes 

tg = 7 X 10" 5 T 3 ^ 2 sec (3-22) 

for temperatures in the eV range. The magnetic diffusion time 
is plotted versus temperature in figure 50. 

3. Density Tailoring in the Nonconducting Regime 
' If t^, is much longer than the characteristic time of the 

til 

process, no thermal conduction occurs. In this regime it is 
relatively easy to change an unfavorable density profile to a 
favorable one. The central part of the plasma is heated, it 
expands and the density drops to a relative minimum on the axis. 
There is no thermal conduction to smooth out the temperature 
profile and consequently the density. Moreover, fairly modest 
energies are required to make a significant density minimum in 
a plasma. 

The simplest estimate of the behavior can be made using the 
work of Raizer 33 on laser focusing and de focusing in atmospheric 
air. This can be applied exactly to weak laser heating of an 
unconfined nonconducting plasma slab so long as the correct 
absorption length is used. Applying Raizer's results, it can 
be shown that for laser pulses longer than the time for an acous 
tic wave to traverse the beam, the density change is given by, 


5 J o /A 

= 1.73 x 10 5 


572 
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max 
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Figure 50. Magnetic diffusion time in a fully ionized plasma 
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where I(r) is the beam intensity profile and 1 is the 
maximum intensity. Here it has been assumed that the absorp- 
tion length, is long compared to the radius of the beam, and 
linearized fluid mechanics is used to evaluate the motion. 
Physically, it is seen that the heating produces a "hollow" 
in the density with the same shape as the beam intensity 
profile, As an example, the density profile that would be 
produced by a laser beam with a parabolic intensity profile 
is shown in figure 51. For practical values, e.g., =» 20 

joules at l.OSp , A = 3 mm 2 , £ = 1/2 , and T e =. 10 2 eV } 
then $p/P 0 =0(1) - This calculation is valuable only as an 
order of magnitude estimate but it does indicate that a 
sizeable density hollow can be produced by a modest laser pulse, 
a. Longitudinal Boring 

Suppose the pre-laser is fired down the length of the 
plasma, then the technique might be called “longitudinal 
boring" since a hollow is being “bored" into the column in 
a longitudinal manner. 

There are basically two ways that a plasma can be bored 
longitudinally. The first way could be used if the plasma 
length is comparable to the absorption length of an available 
pre-laser. The laser beam could be made .very narrow by a 
Couple of lenses (beam shrinking) . In this case the entire 
length of the plasma could be heated simultaneously with a 
very short pulse. The subsequent motion of the plasma 
would produce a density minimum along the axis where most 
of the radiation had- passed. For this method the character- 
istic time is not the pulse length, but the acoustic time— 
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Figure 51. Density produced by weak 
laser heating of a uniform plasma 
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which is the time in which the dynamics will produce a favor- 
able density profile. Then to be in the nonconducting regime 
one must be below the correct t^ = t^ line in figure 4g. 

This corresponds to lower plasma temperatures or higher plasma 
densities. 

Beam shrinking corresponds to a reduction in the parameter 
F j T - (r^/b) 2 . Small values of F allow the beam to be 
trapped for much longer lengths as seen from figures 37 and 
40. of course the smallness of the beam is limited by diffrac- 
tion effects so that there is a limit on the length that a 
beam can be contained by beam shrinking. The diffraction 
limit for a Hd + glass laser is shown in figure 32 where the 
focused laser results of figure 40 are also plotted. For 
example, at n & = 10 18 /cm 3 (C m = 0.001) with b = 1 cm and 
r^/b = .25 , a column of plasma ~ 10 meters long can be irrad- 
iated. The absorption length at 100 eV of 7.2 m is compara- 
ble. So an Nd + glass laser could be used to prepare such a 
plasma for later heating by the main laser. 

R second approach to longitudinal boring uses the favor- 
able effects of plasma motion that arise during longer laser 
pulses. The physical process that will ensue is described 
briefly. When the laser pulse begins, the beam is refracted 
out of the column in a few diameters. But the absorption and 
the higher intensity in the center of the beam tends to heat 
the central portion more strongly. It expands and a very 
weak density hollow is formed right at the end of the column. 
As soon as it is formed, the refraction losses cease at the 
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Figure Diffraction limit for a 1.06^ 

Nd^ glass laser with a focused laser beam 
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hollow, so the beam can travel a few more diameters down the 
column before being refracted out. Stronger heating of the 
central region of the plasma beyond the density hollow occurs. 

It then expands and extends the hollow farther down the column. 

The process continues with the hollow- formation region moving 
on down the plasma column much like a wave. Behind the hollow- 
formation region is a trapped beam in a favorable density pro- 
file. Ahead of the follow- formation region is a rapidly diver- 
ging beam in an unfavorable density profile. The "laser drill" 
as it might be called is shown schematically in figure 53. So 
then a curious coupling of the refraction and the dynamics pro- 
duces a phenomenon whereby the laser actually "drills" out a 
density hollow in the plasma column. The column is then pre- 
pared for the main laser whose beam will propagate down the 

column without refraction loss. 

Xt is clear that the .laser drill will work in the noncon- • 
ducting regime. But it is not certain that it will work when 
thermal conduction effects become important— for they tend to 
prevent the formation of density hollows. The possibility of a 
laser drill in the isothermal regime 'is discussed in section B 4, 

The velocity of the laser drill is an interesting calcula- 
tion to make. This cannot be calculated exactly without solving 
the complete problem in which dynamics and refraction are 
coupled, but a suitable estimate can be made. 

' The length of the hollow- formation region is roughly the 
length in which the beam is refracted out of the unfavorable 

density region, i.e-, Then the velocity of the drill 
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Figure 53. The laser drill 
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must be on the order of Xj. divided by the acoustic time. 
Otherwise the drill would be moving faster than the plasma can 
move to establish a favorable profile-vhich is impossible. 
Hence the governing condition is 


V 


drill' 


X % 7 


2r. 


where a is the acoustic speed in the plasma. So the drxll 
velocity exceeds the acoustic speed by a factor of Xj/ 2r b * 
■This condition is necessary for another reason. If the drill 
velocity does not significantly exceed the acoustic speed, 
then the plasma will squirt out the ends of the column before 
the drill gets moving. 

b. Lateral Boring 

There is another technique whereby a plasma column may be 
selectively heated from the sides to produce a favorable den- 
sity profile. Xf a laser were fired into the column from the 
side, then the denser parts of the plasma would absorb more 
energy. But in an unfavorable profile, the highest density 
is in the center. Hence, if one is in the nonconducting 
' regime, the hotter center region will expand to produce a 
density minimum., a favorable profile. The requirement for 
this regime is simply' t g < <= th , i- *** fche appropriate 
line in figure 49, except now the characteristic distance is 
2b instead of 2r fe . This merely involves identifying the 
lines in' figure 49 with b instead of - A possible 

geometry for lateral boring is shown in figure 54. 
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Figure 5 l K lateral Coring of a plasma column 
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4. Density Tailoring in the Conducting Regime 
The conducting regime is dominated by thermal conduction 
so that temperature nonuniformities are small. This is a majc. 
problem in a plasma with an unfavorable density profile since 
there must be a higher relative temperature near the axis to 

achieve a dsnsity winiitiun\. 

a. Unconfined Plasmas 

The problem can be amply illustrated by solving a simple 
case — weak laser irradiation of an infinite plasma slab in the 
isothermal regime. Assume small perturbations in temperature 
t /T q = 1 + 6* , with 9- « 1 ' where T 0 is the temperature 

of the unheated plasma. Let the dimensionless distance and 
time be R - r/r b , and t = t/t p , respectively. Suppose the 
heating is slow enough that pressure equilibrium is maintained, 
j _ t « t . Then the energy equation is 

i.es., «- a P 

Y+l 501 _ n . I A ( R Mi ) = a- i(R) 

~7~ n R &R V K OR / 


where Y is the ratio of specific heats. The thermal conduc- 

i _ x. /f- is larqe in the isothermal regime 
tion parameter, r\ - ^p'^th ? 

, i „ j. tx. 4 e small for vssk heating 

The absorption parameter a - t p /t h 

The density is assumed to be less than the critical density . 

The dimensionless intensity profile i(R) is the quotient 
of the intensity at a point in the beam and maximum intensity. 
This equation can be solved using Hankel transforms and 


has the solution. 



o» * s* * T 

0'<R,t) = - J IwT (e Y+1 ) sJ o (sR)ds1 , 

{3—23} 

where i{s) is the Hanhel transform of i(R) and J c is the 
Bessel function of order zero- Since the pressure is constant, 
the density perturbation C' “ <P ~ P o >/P 0 is * e;! - at8a to tha 
temperature by =s — 0 1 - Sxamxnatron of (3—23) reveals 
that the initially uniform density profile decays in roughly 
the time ij*t ® 1 (t decay = t^) to a steady state P rofile 
with a hollow the exact shape of the intensity profile. Thus, 

i(R) + 0(^r) (3-24) 

for p'T » 1 . The resulting density profile resembles that 
shown >1. The depth of the density hollow is 

- Ap/p i = Yt _/( y+D t h . 

p 

Thus, it would seem that even in an isothermal plasma, 

density hollows can be produced. But, as soon as the pulse 

/ 

ends, thermal conduction will fill in the hollow in a time of 
order t^ , so one cannot delay firing the main laser into 
the favorable density profile. 

b. Magnetically confined Plasmas 
The conclusion is a bit optimistic when the reality of a 
nonuniform initial density profile is faced. Instead of malting 
a density hollow in a uniform density profile, one must malte a 
density hollow out of what begins as a density "hump". This is 
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a more difficult problem. Complicating the problem is the fact 
that the. magnetic field is nonaiffusive as seen from figure 50. 
The magnetic diffusion time is very slow compared to times of 
interest, so the magnetic lines of flux remain fixed to the 
fluid particles. The disadvantage of this is that the plasma 
has a "memory" , it remenbers its initial density profile and 
is inclined to retain a profile of similar shape. Thus, as 
soon as the laser pulse ends, the nonuniform temperature decays 
away in a time t^ and a density profile resembling the ori- 
ginal one is assumed. 

A discussion of the problem of heating an isothermal plasma 
column contained by a nondiffusive magnetic field is given in 
appendix two. There it is seen that it is impossible to change 
a density hump to a density hollow in an isothermal plasma 
contained by a nondiffusive magnetic field. 

Thus, the question of making a favorable profile out of 
an unfavorable one is simply a matter of making the hollow of 
(3-24) deep enough to over power at least the center of the 
hump in the density profile. An approximate expression for 
the net density profile would be simply the sum of (3-24) and 
the initial profile. This is reasonably accurate under the 
assumptions of (3-24) J It disregards the effect that the over- 
all temperature of the plasma will increase slightly due to 
the heating. 

It becomes even easier to make a density hollow if the 
radius of the laser beam is somewhat less than the radius of 
the column. Then the adverse curvature of the density profile 
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will seem much flatter relative to the narrower beam. This 
effect is demonstrated in figure 55 for a parabolic density 
profile and a laser pulse of 22 joules, a column radius 
b S 1 CP, n Q = 5 X 10 18 /cm- 3 , X 0 - 1.06u, for various beam 
radii. It is seen that the smaller beam radii produce deeper 
but narrower density hollows. The larger beam radii produce 
no density hollows at all. 

Then it is possible to generate a favorable density pro- 
file in a conducting plasma. It is done by continuous energy 
addition to a localized region in the center of the column. 

This is made possible by the fact that when there is continuous 
energy addition, the temperature really isn't spatially uniform 
but may have a density hollow. Therefore there is a possibi- 
lity of having a laser drill in an isothermal plasma column, 
if the absorption length is less than the length of the column. 

It is also possible to use the beam shrinking approach in the 
conducting regime, except the pulse must be long enough for 
•the ensuing dynamics to create a density hollow. 

5. Summary 

There are many varied approaches to achieving a plasma 
column with a favorable density profile. One is to catch the 
plasma during its formation when it may have a favorable profile 
due to collapsing or reflecting shock waves. The other is to 
operate on the guasi-steady plasma column in such a way as to 
produce a. favorable profile. 

The -laser which is used in the operation need not have the 
same wave length as the main laser. In fact the low temperatures 



172 



173 

seem to demand a shorter wave length for the pre-laser , 

' laser irradiation of a plasma column falls generally into 
one of two regimes. The nonconducting regime has no thermal 
conduction and will arise for lower temperatures. The con- 
ducting regime is dominated by thermal conduction and occurs 
for higher temperatures, 

Xf the plasma column is about as long as the absorption 
length of the pre-laser, then the whole column can be bored out 
in the center by beam shrinking, mating the laser beam narrow 
enough to travel the entire length without being refracted out. 
If beam-shrinking doesn't work, then the laser drill can be 
used whereby the shape of the intensity profile actually drills 
out a density 'hollow along the axis . The plasma column can 
also bo bored by heating from the side (in the nonconducting 
regime) . 

There are certain peculiarities of the conducting regime. 
It is' possible to generate a density hollow but it persists 
only as long as the laser pulse. After the pulse ends, the 
"memory" of the nondiffusive magnetic field returns the plasma 
to an unfavorable profile much like it had in the beginning- 
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C. APPENDICES 


Appendix 1 

Behavior Near the Axis and Near Caustics 
Applying directly the results of geometrical optics leads 
to infinite electric fields at both caustics and at the axis. 
They arise at caustics since the rays envelope. They arise 
at the axis simply due to the cylindrical symmetry, a wave 
collapsing to the axis will generate a singularity in the 
intensity {at least in the linear case) . But singularities 
in electric fields do not occur in nature and thus there must 
be an effect that limits the magnitude of the electric field. 
The failure of geometrical optics can reasonably be expected- 
one cannot apply it to phenomena occuring in lengths shorter 
than a wave length. 

The study of these singularities begins with an examina- 
tion of the full linear equations. Eliminating the magnetic 
field intensity from Maxwell's equations, an equation for the 
electric field can be found. 


0 » 


Vx(7xE) + p 


bj 

bt 


+ 


1 & 2 E 



(3-25) 


*4 ^ 

where E is the electric field intensity and J is the current 

density. Consider only the current due to electrons, J = " en e v e 

where v is the mean electron velocity which is governed by 
e 

an electron momentum equation. 



175 


(3-26) 


m.. 


bv 
e 

bt 


-> 

-eE 


-5 

m v 


e e 


er 


The last term in (3-26) is the electron— ion collision term. 

It assumes the electrons lose their ordered oscillatory motion 
to random motion due to collisions in a time T ei * ' T ei ls 
the electron— ion collision time. Considering only oscillatory 
fields, E = E{r,z)e -:LtUt , then.all time derivatives can be 
replaced by the factor -iiu . Then (3-25) and (3-26) combine . 
to yield 

-V 2 E - V(V-E) + Tc 2 E = 0 , (3-27) 


where k is the magnitude of the wave vector 


k 2 (r) 



2 , , 2 
a) +lA e 


(l-i/ujT)1 


which depends only on r . 

If attention is confined to the case where the electric 
field is in the azimuthal direction, then (3-27) becomes 


b 2 E . 1. JL 
,_2 + r 6r 



+ k 2 E = 0 


(3-28) 


Suppose now that in addition to k being independent of z , 
that E is also independent of z except for periodic spatial 
oscillations. This assumption is not exactly true in a real 
plasma column irradiated from the end by a laser. In that case, 
rays entering the column at different points will trace some- 
what different shaped paths. This model has all rays tracing 
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the same paths — so will only be locally valid. In other words, 
the behavior of a particular ray in the real case is being simu- 
lated by a solution in which all rays trace, the same shaped 
path as that particular ray. 

Then one can take E = F(r) exp[±k(r) sin0{r) z] where 
6 (r) is the angle between the wave vector and the r axis. 

Then (3-28) becomes, ' eliminating the dependence on z s 

i 4 - (r 4^- ) + [k 2 (r) - k 2 sin 2 9 (r) ]F = 0 , 

r dr v dr o 

and 

k(r) sinB (r) = k 0 sin6 o , 

where k Q is the wave vector at r = 0 and 0 O is its angle 
with the r axis at r = 0 . (3-30) is the condition that 

must be satisfied to eliminate the dependence of F on z . 
(3-30) also gives the angle of the rays at a particular r 
and gives the ‘‘turning point" where a ray is reflected. Figure 
56 plots the radial position of the turning point versus the 
angle of the ray as it emanates from the axis for the case of 
the parabolic density profile (3-14) with C m = O- 1 * 

The turning point in this approximate model corresponds to 
the point where a ray borders the caustic in the real problem. 

In the real problem the ray touches the caustic a short distance 
after achieving the maximum distance from the axis not right 
at the maximum point. 

Applying (3-29) to a parabolic density profile (3-14) 


(3-29) 

(3-30) 


yields. 



17? 
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2 

d tr g) + * 2 4 {1 ' C J CcOs29 0 “ ^ + ° ( ^ ei 

G o 


•)3p = 0 . (3-31) 


dr 


lt is seen that the coefficient of F in (3-31) becomes small 
at two places? at ’ r = 0 it actually . vanishes, and near r = b 
cos 0 o ’ (the turning point) where the coefficient becomes small 
(of order 1 /ujt) - Hear the turning point, r = b cos 9 o , the 
solutions are Bessel functions of one third order, 


F « const y 1 ^ (Jjy 3 (y) + J -l/3^ ^ ’ 


with 

3/2 

2 (1-C ) cos9 r r, COS o i_~j 

„ _ bm m — :: % — I (cos 8 - r-) o . 2 Q uutJ. 

y *- c 2,. in 2 fl u 0 b cos^e +C m Sin 0 

0 COS 0 +C sin a 0 0 m 0 


where b cos 8 q - r « b 


This type°of solution has been discussed by Ginzburg. 21 
Physicaliy it represents a standing wave pattern near the turn- 
i»g point ana to the superposition o£ the incoming and reflected 
rays. There is not a real singularity since y J -l/3 w 
approaches a constant value as y - 0 . The electric field 
intensity near the turning point will exceed that™ nonsingu- 
lar regions by. a factor on the order of <b*/c„) 

, • _ _ n f3-31) approaches Bessel's equation 

Near the axis, r - 0 , ^ 

of zero order so that the solution will be for r « b , 


F * const J (f* O + const Y 0 (? 
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where J and Y o are the Bessel functions of zero order of 
the first and second kinds respectively. It is Y q that has 
an actual singularity. Geometrical optics predicts a 1/r 
singularity in the intensity or a l/r^ singularity in the 
electric field. For large values of the argument of Y q , 


1/2 


(s ./gr) . (-Jfr-) 

' c - m ' rrujr,yi-C 0 


5rtN 
" 12 / 


which (aside from the oscillatory part) varies like 1/r . 

and agrees with geometrical optics. But geometrical optics 
fails for distances smaller than a wave length, r < c/w , 
i.e., when the argument of Y 0 becomes small. Then, 


Y o 




which is a much weaker singularity— though nonetheless still 
present. Physically this means that most of the wave is 
reflected within a wave length of the axis and only a small 
part ever converges all the way to the axis to generate a 
singularity. 

Of course the axial singularity is not a real one. Non- 
linear effects will stop the singularity at some point. More- 
over the significance of the singularity is doubtful anyway. 
Only for small distances much less than a wave length) will 
the electric field be very large. But the electron mean free 
path (for collisions with ions) is never much less than the 
laser wavelength for thermonuclear plasmas. 


There can be no 
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heating without electron-ion collisions, hence the temperature 

couldn't have a singularity even if the electric field did. 

There is another effect that has been observed along the 

axis of a laser beam. That is the collapse of part of the 

beam into an intense filament. 34 This has been attributed to 

3 5 36 

nonlinear effects of intensity on the index of refraction. ' 
This effect — which is not well understood may hold possibi- 
lities for the trapping of very narrow laser beams. 

Appendix 2 

Heating of Conducting Plasmas 
Contained by Nondiffusive Magnetic Fields 
A solution has been found for the heatrng and expansion of 
a plasma column which is in the conducting regime. The heating 
is slow enough that pressure equilibrium is maintained, 
t « ,and the heating is conducting, t fch « t & . Thrs 

' requires a plasma temperature above the . line t^ = t a in 
figure 47 and above the line t Q = t^ in figure 49. 

In a cylindrical magnetic field, the magnetic field, B , 
is related to the magnetic flux, $ B , by 

D _ 1 __B 
. B “ 2trr 5r 

Since the field is nondiffusive, is fixed to the fluid 

particles and moves with them. The displacement function for 
a particle is r = X(r„,t) , where X is the location at time 
t of the particle that was at r o when t = 0 . Then § B 
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is a function of r Q alone, so that 


B = 


a$ B /dr o 

„ £>X 

x '5r~ 

0 


Than the magnetic pressure p„ - E /2p becomes 


P B P B 0 ° 


2/x -£2L 


(3-32) 


where p (r } is the magnetic pressure profile before the 
heating; t = 0 ♦ 


Nondimensionalizing the variables, the pressure scale is . 
the total pressure before the heating. The temperature scale 
is the initial temperature which is spatially uniform like 
the pressure. The density scale is the maximum initial den- 
sity which is on the axis. Then the equation of state is the 
pressure, fl , equals the sum of the magnetic pressure II B 
and the thermodynamic pressure C 6 » 


H = C8 + n B 


(3-33) 


Initially, the equation of state is, 

1 = c 0 +n Bo 


(3-34) 


The Lagrangian continuity equation relates the density to its 
initial value using the displacement function, 


where 

Then, 


i (r ) 

o o 

combining 



o 


is the initial 
(3-32, 33 ; 34, 


=o o 

density profile (at 
35) , an equation for 


t 


(3-35) 

= 0) . 

X results. 
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6r 


£ (x - (l-£ ) r 2 = 0 

o^o\ br„ / o o 


(3-36) 


Since n and 0 are spatially uniform and depend only on 
time, (3-36) can be integrated exactly and has the solution 


X 


2 



+ kl + |(1”C 0 )3 t 0 dr 0 


(3-37) 


The density profile £ can be written using (3-35) and (3-36) , 


c 


2E 

T* 


1 + 


/ 


4n 1 ~ g o 

8 2 C 2 


(3-38) 


The pressure, II , and temperature, 8 , are functions of 

time which depend on the laser power and on the character of 
the confining magnetic field. 

Unfortunately, (3-38) gives £ as a function of r 0 , 
but one is interested in £ ( r ) ♦ it is then necessary to 
invert the displacement function found in (3-37) ; t “ X (r,t), 
which probably cannot be found explicitly in most 

examples of interest. But at least the solution is available 
in quadrature form. 

The solution has been calculated for the case of a para- 
bolic initial density profile, £ o = 1 - (r Q /b) . The dis- 
placement function using dimensionless variables, y = r/b = 
X/b and x = r /b , is 
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(3-3! 


where A » ^ — 2 . 

Tin approximate expression for the density profile near 
the axis is 

C =|(l -~-y 2 + 0(y 4 )) (3-40) 

It is seen that it is impossible in an isothermal plasma to 
change an unfavorable density profile to a favorable one. 

(3-40) shows that it is still unfavorable with maximum C 
on the axis. The only way then to get a favorable profile 
•is-- to have continuous localized laser heating as described 
in section B 4. Then a small density hollow can be maintained 
as long as the laser pulse lasts. 



Chapter Four 
DIFFUSION WAVE PHENOMENA 


A, THE BLEACHING WAVE 

One of the really interesting results of chapter II on 
laser heating of a stationary plasma was the appearance of 
the so called bleaching wave. This phenomenon differs signi- 
ficantly from the usual linear absorption. In the linear 
absorption problem, light is attenuated by a substance whose 
absorption coefficient is a constant. The result is that the 
intensity decays in simple exponential fashion. 

In the case of laser plasma heating, the absorption 
coefficient is not a constant but decreases with increasing 
electron temperature (K v ~T e ~ 3/2 ) . Generally, K v is quite 

large in the unheated plasma so that tire radiation is atten- 
uated in a very short distance, i.e„, the plasma is opaque 
to the radiation. If the absorption coefficient is large, 
then all the laser energy is added to a very thin layer of 
plasma. The energy addition causes the temperature to increase 
rapidly, so that K y decreases— even more rapidly. Thus, the 
layer becomes successively more transparent as time proceeds. 

Then the overall behavior is characterized by the propa- 
gation of a bleaching wave into the plasma. Successive layers 
absorb the radiation and become transparent. In essence, 
the radiation bores its way into the plasma. Behind the 
wave, the plasma is highly excited and nearly transparent, 
so that the radiation is only weakly attenuated. In the front 
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of the wave is a region of rapid attenuation of the radiation 
and consequently rapid "heating". Ahead of the wave is a 
weak exponentially damped beam in a "cold" plasma, closely 
resembling a linear absorption problem. 

The term bleaching wave is an appropriate name for this 
phenomenon. As the wave moves, it bleaches the opaque plasma, 
making it transparent to the laser radiation. 

It seems intuitive that any material having an absorp- 
tion coefficient that decreases under continuous radiation — 
will exhibit a bleaching wave behavior. This is indeed true 
and several illustrative examples are discussed. However, 
first a presentation of a general example is given to demon- 
strate the bleaching wave. 

The radiative transfer equation for radiation propaga- 
ting in the tx-direction and neglecting thermal radiation is 


bI v 

0 = ITx + Vv 


(4-1) 


where I v is the radiation intensity. The absorption coeffi- 
cient is some function of the temperature and the density, p , 
K _ K v (p,x) . The energy equation for a stationary substance 


is 


bT 


pW c v I = - 


bi^ 

bx 


(4-2) 


where the right member represents the divergence of the energy 
flux (intensity), C y is the constant vo'lume specific heat, 
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and p(x) is the density. The initial condition is the 
spatial temperature distribution at time t = 0 ; T(x,0) - 
T (x) . The boundary condition is the intensity entering 
the substance at its edge (x = 0); X v (0 J t) - (t) which is 

the intensity delivered by the laser. 

Combining (4-1) and (4-2) yields a second order partial 
differential equation 

o - & [ t & !pT) 3 * & (pT1 <4 - 3 ’ 

In general, the absorption coefficient is a function of 
p and T : K y = K y (p,T) . This can also be expressed as a 

function of p and pT : K y = K*(p,pT) . Over any given 

range of the variable pT , the anti-derivative of 1/ K V 
with respect to pT exists such that 

61, ( p . pT) 1 (4-4) 

6(pT) K*(p,pT) 

X, being the anti-derivative. The anti-derivative is not 
unique because if any function of p is added, (4-4) is 
still satisfied. Then one can say that L(p,pT) = L Q (p,pT) 

+ F(p) , where L is chosen as the principle anti-deriva- 
tive and Ftp) is an arbitrary function. 

Based on the existence of the anti-derivative, and the 
independence of p on time, (4-3) can be rearranged to give. 
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0 - H [ TS + F > + pT ] (4-s; 

Using the initial condition T{x,Q) * T^ (x) , (4-5) can be 
integrated with respect to time to yield. 


0 ** h [ L o (p > pT) - L o <^P T 0 >] + P T - pT o (4 “ 6 > 


Notice that the arbitrary function P(p) has disappeared in 
the application of the initial conditions. 

The boundary condition can be written as 


ft) = ^ M j 


which. in turn can be integrated, applying the initial condition 
to yield 

t 

[i, o (P,PT 5 ~ Vp,pt )] i X=G j V*'^' * 

( 4 - 7 ) 

The first order partial differential equation (4-6) together 
with the boundary condition (4-7) comprise a well posed pro- 
blem. These are the general equations governing temperature 
in a one-dimensional medium irradiated by a laser. 

If the absorption happens to be linear (K y independent 
Of pT), then is just pT/K v • Then the familiar 

exponential decay temperature profile will appear. For 



example, if the initial temperature is zero, density, laser 
intensity, and absorption coefficient are constant, then the 
temperature profiles are 




The temperature is exponentially damped as one moves into the 
plasma, and rises with a uniform time factor as time moves on. 
This illustrates the basic character of linear absorption. 

There are two types of nonlinear absorption; K y decreas 
ing or increasing with temperature- The former rarely occurs 
in nature and does- not produce a bleaching wave except in 
the case of a focused beam. The latter leads to the bleach- 
ing wave. 

' (4-6) cannot be integrated in general. Three physical 

examples are presented for which solutions to (4-6) have been 
found, baser heating of a stationary plasma was presented in 
chapter II . Laser irradiation of a resonant two level sub- 
stance, and the laser induced breakdown wave are presented 
in this section. 


, 1. Laser Irradiation of a Two Level Substance 
Generally an atom or molecule has many discrete energy 
levels in its internal structure. A molecule has one. or more 
vibrational modes, each with its own set of energy levels and 
the corresponding quantum numbers. An atom also has a com- 
plex set of energy levels and corresponding quantum numbers. 
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Transitions by individual particles from one energy level to 
another are made by several means. For example , a particle 
in a lower energy state of energy E ^ can be raised to a 
higher energy state E^ by the absorption of a photon of 
of exactly the right energy; hv «* B 2~ E 1 * Or, a particle 
in the higher' energy state when bombarded by the same photon, 
will drop to the lower energy state, emitting a photon of the 
same energy — so that two identical photons result, which is 
called stimulated emission. ' 

Consider a substance that has more particles in the lower 
energy state that is being irradiated by a laser of exactly 
the right wave length so that h v = E 2~ S 1 * ^en, on 
statistical average more photons are being absorbed than are 
being emitted by stimulated emission. This is simply because' 
more particles are in the "absorbing'* state than are in the 
'’emitting" state. But, as the beam is absorbed, more and 
more of the particles are raised to the upper level, so that 
the number in the upper level approaches the number in the _ 
lower level. Now, the absorption is much weaker since the 
ratio of absorbed to emitted photons is reduced to near one. 
Eventually, the numbers in the two levels will be equal so 
that there is no absorption at all, because the number being 
emitted is equal to the number being absorbed. 

It is immediately seen that this is exactly the same 
phenomenon that occurs in the laser heating of a plasma. The 
local absorption coefficient diminished with time under 
continous irradiation. 
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This sort of phenomena is known as saturable absorption, 
a type of nonlinear absorption. It has been studied exten- 
sively since it is this sort of substance that makes a laser 
or a laser amplifier. The only difference is that for a 
laser, there are more particles in the upper energy state 
than the lower so that a beam of .the right wave length is 
amplified rather than being attenuated. 

The solution to this problem was first worked out inde- 

37 38 

pendently by Bellman, et. al., by Frantz and Nodvik, and 
39 

by Schultz -Dubois. Each worked out the general time depen- 

dent solution in a two level substance. These papers neglected 
spontaneous emission and thermal relaxation from the upper 
energy level which is appropriate for short time scales. 
Bellman, et. al., pointed out the appearance of a phenomena 
whereby the radiation "bores" its way through the substance 
with a certain velocity which is simply a bleaching wave. 

Later, Siegmann 40 worked out the solution neglecting the 
transient term in the radiative transfer equation. He obtained 
the same results as in the original studies but in a much 
easier manner with more lucid results according to his claim. 

Presented below is a simple solution of laser irradiation 
of a two level substance. This calculation is similar to 
the approach of Siegmann but it is instructive in that it 
clearly demonstrates the bleaching wave. 

The radiative transfer equation has the usual form (4-1) . 
For this problem the energy equation takes a slightly different 
form, 
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hv 




{4-8) 


where n is the number density of particles in the upper 
state and is analogous to the temperature in {4-2) . h is 
Planch's constant and v is the frequency of the radiation. 
Of course the sum of the number densities in the upper (n u ) 
and lower (n^) states in a constant (n) . 

In a two- level substance, the absorption coefficient is 

K v = a ^ n r n u ) (4 " 9) 

where a y is the absorption cross section of the radiation. 
The initial condition is 


n u (x,°) - n uo {x) 


(4-10) 


The boundary condition is l v (0,t) = (t) 

combine to yield. 


(4-8) and (4-9) 


0 hx \ t>t ) + 6t 


(4-11) 


1/K^ has as its principle anti-derivative, - (l/2a v ) log (n-2n y ) 
so that (4-11) can be integrated with time. Performing the 
integration and applying the initial condition yields. 


0 =*~-or— *s~ [log (n-2n ) -log (n-2n )] * n -n 


2o v 6x 


U Uo 


(4-12) 
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This is in the same form as (4-6) with 
The boundary condition is 

- [log(n- 2 n u )-log(n- 2 n uc> ) 3l x=0 = 

which is analogous to (4-7) . 

A particularly simple solution is available for the case 

n = const. Solving (4-12) for this case gives the solution, 
uo 

n = T1 0 + Js(l-n o ){l+tanh[(i-ri o ) (V §) ^ » (4-14) 

where g " 2 n u /n is the dimensionless measure of n Q . 

5 _ x/ix f and t = t/At are dimensionless distance and time 
respectively. The length scale is Ax = 2 / 0 v n an<3 the t:Lme 
scale is At = hv/2a v l o , I 0 being the average incident 
intensity. The position of the wave as a function of time 
is given by 

t 

where I = I /l . The number density is seen to be a 
O Vo O 

hyperbolic tangent function. It raises g from ti 0 ahead 
of the bleaching wave (f»? w (t) ) to 1 behind the wave 
(5 « § (t) ) . When n = 1 , the substance is said to be 

"saturated" . 


pT replaced by n u . 


f I (t')dt' (4-13) 

V Vo 
0 
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The expression for ? w (t) is composed of a steady state 
term and a transient term. Considering only the steady state 
term, the velocity of the bleaching wave is given by 


v 


I (T) 
o 


Ax 


w 2(i-n 0 ) At 


The velocity of the bleaching wave can also be found by using 
the conservation of energy. This is done by requiring that the 
energy flux delivered by the laser (energy per unit time per 
unit area) must be equal to the energy required per unit 
volume to bleach the substance times the velocity of the 
bleaching wave. The calculation is made possible by the fact 
that essentially all of the laser energy is added to a thin 
region, the “front 1 ' of the bleaching wave. However, the con- 
servation of energy cannot be used to calculate the velocity 
of the bleaching wave in a plasma. Fo,r this case the energy 
continues to be added to the plasma far behind the actual 
front of the bleaching wave. A plasma never saturates com- 
pletely as does a two level substance since an infinite plasma 
temperature is required for perfect transparency. 

The thickness of the bleaching wave is roughly the dis- 
tance in which the argument of the hyperbolic tangent makes 
a change of one, so that 




tv 


1 ~ T >o 


Poe the particular example of n 0 = 0 (no molecules initially 
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in the upper state) and a uniform laser- pulse, Z 0 = 1 , the 
solution becomes 

r\ = %Cl + tanh(S w -§)3 , 

§ w = hr + fn(l-e _T ) 

This example is shown graphically in figure 57 . Also shown 
in the same figure are temperature profiles for the heating 
of a uniform stationary plasma by a uniform laser pulse. While 
certain features differ between the two phenomena, a bleach- 
ing wave is seen to appear in both. Ahead of the wave there 
is apparently little radiation since the variable (g or 9 e > 
changes slowly. In the front of the wave, the variable rises 
rapidly, due to strong irradiation and large absorption 
coefficient. Behind the wave, the variable rises much more 
slowly due to bleaching as the absorption coefficient goes 
to zero. 

2. Baser Induced Breakdown Wave 
For many conditions, an un— ionized gas such as air is 
essentially transparent to laser radiation. At least the , 
absorption lengths are quite long compared to absorption 
lengths that arise in partially ana fully ionized gases, or 
from resonant line absorption (section A 1) . Nevertheless, 
an extremely intense laser beam can cause breakdown of the 
gas. This phenomenon was first observed experimentally by 



195 


laser Irradiation of a two level substance 

number density 

in unper state = 0 

*1 
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Damon and Tomlinson and by Meyerand and Haught. 

Studies of the mechanism of laser induced breakdown show 
that the breakdown begins with either stray free electrons in 
the gas, or with electrons produced by the multiphoton effect. 
The multiphoton effect is due to the simultaneous bombardment 
of an atom by several quanta of light, with enough total energy 
to free an electron. Once an electron is free, the large 
oscillating electric field causes the electron to collide with 
atoms and free more electrons. The cascading process pro- 
duces a large number of electrons in a very short time. A 
clearly defined threshold intensity characterizes this break- 
down. Several theoretical investigations of the breakdown 

43 44 45 46 47 

phenomena followed the experimental discovery. ' ’ ’ ' 

The dynamical behavior following the initial breakdown 

was' seen to exhibit unusual effects. The expansion of the 

ionized region was not uniform in all directions. Rather, 

it moves much more rapidly in the cone of the focused laser 

beam. This phenomena was first explained by Ramsden and 
.48 

Savie and was shown to be a radiative detonation wave, that 
is a detonation wave whose fuel is the energy absorbed from 
the laser beam. 

If the intensity of the laser beam appreciably exceeds 
the threshold intensity for breakdown at the focus, then 
another phenomena can occur. In this case the intensity will 
also exceed the threshold over a certain length of the beam 
and breakdown will occur at these points also. The gas in 
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the breakdown region is stationary and the wave like propaga- 
tion of the region towards the laser source is due to the 
breakdown limit being reached at successively longer times 
for points closer to the source. The speed of this "break- 
down wave" may exceed the speed of all of the possible gas 
dynamical processes. The breakdown wave was studied by 
Raizer 49 and its speed was calculated by a phenomenological 

model. 

It is shown here that the breakdown wave is governed by 
equations which give the bleaching wave solution. The only 
variation is due to the non-one-diraensional geometry which 
causes the wave to move backwards, i.e., moving towards the 
laser source as contrasted to the usual sense of a bleaching 
wave moving away from the source. Otherwise one sees all 
the' usual features of a bleaching wave. On one side of the 
wave (behind) is an opaque plasma. On the other side (ahead) 
is a transparent un-ionized gas. The heating of the gas is 
done by absorption of the laser radiation— most of which 
occurs very rapidly at the front of the wave. 

The nature of the breakdown wave as a bleaching wave is 
seen by considering a solution to a simple breakdown model. 
The radiative transfer equation is given in terras of the 
power , W , 

SB. _ KW = o < 4 - 15 > 

The opposite sign from the usual equation arises since the 
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beam is assumed to travel in the -x direction as shown in 
figure 58. Figure 58 also shows the geometry of the focal 
region, fit any point the intensity is related to the power 
by i s w /h (x) , where A (x) is the cross sectional area 

of the beam. 

The energy equation is 

oMxlc M = M (4-16 

. * 

where pfi (x) is analogous to the density in (4-2) and the 
power W is analogous to the intensity. Combining (4-15) 
and (4-16) yields. 



which is clearly analogous to (4-3) . 

This problem can be easily solved with the following two 
assumptions. First the absorption coefficient is given by 

K = } H(n e -n et ) (4-17) 

where & is a short absorption length, n efc is the electron 
number density for breakdown, and H( ) is the Heaviside 
function. Thus there is no absorption for electron densities 
below the threshold, and strong absorption above the threshold. 



focal area 


Figure 58. Geometry of the focal region 
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Following Raizer, 49 the second assumption is that the 

V 

electron number density grows at a rate proportional to the 
intensity times the number density. 





t*W. 


max 


) 


n W 


(4-18) 


where A f is the focal area, V? max is the maximum power 
and t* is the characteristic time of breakdown at the focal 
area. 

Since the absorption coefficient depends on the electron 
density, it is easier to solve the problem in terras of n Q 
and then find the temperature. 

Solving (4-18) for W and applying to (4-15) to eliminate 
W yields 


0 .iif 

u ~ K bx x n Ot / n 


bt 


(4-19) 


Applying the absorption coefficient (4-17) to (4-19) 
yields 


0 = it lo * n e> " T H(n e~ n et> £ (A(x) 109 "e? 


■ The initial conditions are 


(4-20) 


n Q (x,0) - n £ 

T(x,0) 


(4-21) 
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where n fio is a small residual electron density due to chance 
free electrons in the gas, or due to electrons produced by 
the multiphoton effect. 49 The boundary condition is 

W(“,fc) = W L (t) , (4-22) 

where W T <t) is the' time history of the laser power. (4-20), 
(4-21) , (4-22) , plus (4-18) relating n e and W form a well 
posed problem. (4-20) can be solved in two parts; for 
n < n et and for n e > n et , and these two parts can be 
matched. First solve for n Q < n efc so that the Heaviside 
function is zero. Then (4-20) becomes 6 2 (A (x)n e )/bxDt = 0 , 
for which the solution is 

A(x) log n e = a x (t) + b^x) , (4-23) 

where a ] _ and b 1 are arbitrary functions in t and x 
respectively. Applying the initial condition (4-21) yields 
0 = 3 ^( 0 ) + • But fchis can onl y be true i£ is 

a constant. Without loss of generality, this constant can 

be taken as zero so that b x (x) = 0 = a i(°) • 

Applying the boundary condition requires a little intui- 
tive thinking. Consider that for n g < n 0t there is no 
absorption. Consider also that the successive points of 
breakdown (where -n e - n 0t ) will have to move in the +x 
direction as time proceeds. Then, if at some point in time, 
breakdown has not occurred at a point x , tnen it surely 
has not occurred at any place between x and the laser. 
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Then the intensity at x must be the undiminished intensity 
delivered by the laser. . So then the boundary condition (4-22) 
actually becomes 

W(x,t) = W L (t) , for n e < n et . 

Applying this to (4-18) using the solution (4-23) yields 


da 


dt 


r = ( 


n et A f 
^max 


) V« 


which integrated using a^fO) - 0 becomes 

t 


= (^r~) I V t,)dt * • . (4 " 24) 


max 


Then the electron number density is 


n. 


n 


00 


A • 

■“fiHT-ffiff V t '> at '5 W M- 25 > 

max 


Since the absorption coefficient is zero for n 0 < n gfc , the 
temperature will not rise at all, but remains at T Q in this 
region. 

The location of the breakdown wave x w (t) is at the 
point where n & = n gt which using (4-25) is given implicitly 
by 


_L w 

A- max 
f 


t* 



(2-26) 
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'Of course this is precisely the point where the solution 
(4-25) becomes valid. 

Now the solution for n & > n gfc must be constructed. 

Then (4-20) becomes 

0 = shz (a(x) log - A ft A(x) log n J 

This equation can immediately be integrated with respect 
to time, resulting in. a first order linear differential equa- 
tion in x for (A(x) log n J . Solving this equation 

/ 

■ yields 

I 

log n e » a 2 {t) e x/ * + b 2 (x) (4-27) 


^ At this point it is possible to go directly to an equa- 
tion for the temperature behind the wave rather than solving 
for n e . Using (4-18) and (4-27) , the laser power at a 
point behind the wave is 


W „ ( t - %— ) a| (t) e x/l 


The boundary condition is W(x w (t),t) = W L (t) so that 


4<t> = (WH e 


-x^t) l 


max' 
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and 

W *= W L (t) exp[j (x-x w (t) )] (4-28) 


Then the temperature can be found by applying (4-28) to the 
energy equation (4-16) and integrating with respect to time, 

T = -pcW f V' ,Mt ' + “ nSt 

r V o 

The constant is found by applying the condition that 

T(x (t),t) = T q . Evaluating the constant, the equation for 

the temperature can be written, 

~x w (f) ) ]dt* 
(4-29) 

where x (t) is given implicitly by (4-26) . 
w 

Another form of (4-29) which is more amenable to calcu- 
lation involves the inverse of the function , 



T " T o + piC A(x) 




W T (t') exp[-| (x-x w (t*) )]df 


(4-30) 


where 


V*) 


is given implicitly by 

M 


n. 


i * pf 

-i- w t* log 
A f 'max n eo 




A(x) 


W T (t‘) <it‘ 

Li 


(4-31) 


The wave motion is clearly seen in (4-29) where the Heaviside 
function in the integral implies there is no heating until the 
wave arrives , at which time the temperature rises . 
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The velocity of the breakdown wave can be found by taking 
'the time derivative of (4-26), 


vi 

, °et “max 
log - — 

X=JC w (t) eo 

This velocity could also be calculated by applying conservation 
of energy. 

The wave is more clearly seen by calculating an example. 

Consider a uniform laser pulse so that W = W H{t) , T = 0 , 

2 

and a focal region where A{x) = + frj?;) > i* e *> a conical 

beam shape, where r - is the radius of the focal area. For 
this case the solution is 



T = 


W et* 

max 


x' -X) 

(1 (tj)e * 

r f 


U 


'[1 - e " i(Xw<t) " X> ]3 


These temperature profiles are plotted in figure 59. for various 
times. The temperature profiles have the same wave character 
as the laser heated plasma and the laser irradiation of a two 
level substance. The only major difference is that the bleach- 
ing wave moves backwards. This is due to the geometry — the 
focusing effect on the laser intensity. 
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-Figure 59. Breakdown wave 
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Thus three examples of the bleaching wave have been pre- 
sented which bear similar characteristics. In the first two, 
laser irradiation of a layer of substance in time produced 
a significant decrease in the absorption coefficient. The 
result was a boring effect of the laser beam into the sub- 
stance. In the third example, laser irradiation produced in 
time a sharp increase in the absorption coefficient. This, 
together with the focused geometry caused a reversed boring 
effect to appear. The penetration of the beam diminished 
with time’. All three of these phenomena show similar pro- 
files in the dependent variable, as seen in figures 57 and 
59 . In the next section, another phenomenon will demonstrate 
this same wave like behavior. 
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B, THE THERMAL WAVE 

There is another class of phenomena which seems to have 
a different physical base than the bleaching wave but whose 
solutions exhibit a vary similar behavior. Thare are several 
manifestations of this in physics as will be discussed in 
section D, but in this section, the discussion will be con- 
fined to one example, thermal waves in media with nonlinear 
thermal conduction. 

In thermal conduction problems, one usually assumes a 
constant coefficient of thermal conductivity. This is moti- 
vated by a desire to make the equations tractable-~linear in 
this case so that superposition and all the other techniques 
associated with linear equations can be used. Yet, there is 
no substance found in nature which has a constant thermal , 
conduction coefficient: it always depends on the temperature. 

Only for small changes in temperature will linear analysis 
apply— and then only approximately. 

Most real diffusion problems are nonlinear and numerous 
examples have been cited by Cole . Some of these examples 
will be mentioned in section D. A thorough review of the 

51 

mathematics of nonlinear diffusion has been given by cranX. , 

The numerical approach to nonlinear diffusion problems has 

52 

been discussed by Richtmeyer. 

A typical nonlinear thermal conduction problem is governed 
by an equation of the form, 

»=v It - to ( K(T1 ) - 0 


{4—32) 
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where T is temperature, K(T) is the coefficient of thermal 
conduction, and Q is a source term if needed. Suppose that 
at the edge of a substance, energy is being added. Then ther- 
mal conduction will transport the energy from the edge to i deeper 
layers of the substance. At any point in time, the temperature 
profile will be monotonically decreasing with increasing dis- 
tance into the substance. If the thermal conduction coeffi- 
cient were constant (linear diffusion) then any temperature 
profile would be an exponentially damped curve. 

However, if the thermal conduction coefficient is greater 
for higher temperature, then heat will conduct more easily 
in hotter regions, i.e., nearer the edge. Thus the hotter 
regions will require smaller temperature gradients to conduct 
a certain amount of heat than cooler regions. Hence the 
temperature profiles will exhibit the humped shape seen in 
the laser heated plasma of figure 57. 

This is clearly seen if a particular example is worked 
out, that of a fully ionized plasma. In a fully ionized 
plasma, the thermal conduction coefficient is proportional 
to T^ 2 , 53 Suppose the source term in (4-32) corresponds 
to a constant heat addition rate at the edge of the plasma, 

Q = Q 0 6 (x) R(t) , where 5{x) is the delta function. If 
the initial condition is T(x,0) = 0 , and the boundary con- 
dition is T(”,t) = 0 , then (4-32) can be written as an 


ordinary differential equation using similarity methods. 

~ ~7 /g ~ 

The similarity variable is r = x / t > where x and 


t 
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are suitably chosen dimensionless distance and time, respec- 

ai /v2 /o 

tively. The dimensionless temperature is T = t ' f (r|) . 
f (n) is governed by 



<v 

where K ‘ is the dimensionless parameter representing the 

thermal conductivity. The boundary conditions are f{®) = 0 
7/2 

and £i| (o) = - l . Then for length and time scales for 

dr] 

which the conductivity parameter is large, K » 1 , . an approx- 
imate solution can be calculated, 


£(n> 



This solution becomes invalid for ln -T ? 0 l <<: 1 » For h > tig 

the solution is £ = 0 . Near q = t] 0 > a boundary layer 

approach is necessary to give the details of the solution. 

•n is a constant which can be evaluated by integrating over 
o 

x to get the total energy added, and equating this to the 
time integral of the source term. Conducting the integrations 
yields n 0 - (9/7) 7 ^ . Then the temperature profiles are 

T = (T) 0 t ' ~ >0 


Several of these profiles are plotted in figure SO for a 
fully ionized plasma. The "boundary layer" at the front of 
the wave is not shown. 
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Figure 60. Nonlinear thermal conduction 
in a uniform plasma slab ■ 
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The presence of wave like behavior is clearly seen. This 
has been appropriately called a thermal wave by Zel'dovich 
ana Rai 2 er. 54 The resemblance is even more striking due to 
the differing character of the governing equations. The thermal 
wave is governed by an equation of the form of (4—32) whrch 
is parabolic. The bleaching wave is governed by an equation 
Of the form of (4-3) which is hyperbolic. This resemblance 
is carried over to the linear case where laser irradiation 
ana thermal conduction both produce exponential decays in the 
temperatures produced. 

The resemblance of the thermal wave to the bleaching wave 
is not presented as a mere curiosity. Nor is it presented to 
introduce a coincidental resemblance between two physical 
phenomena. Rather, the goal is to demonstrate that these two 
- phenomena — which arise out of ’widely differing equations-- 
actually are due to two opposite limits of a physical process, 
i.e., they both fit into a general class of "diffusion waves" . 
There are many physical phenomena which exhibt these diffu- 
sion waves. This contention is presented and supported in 


the next section. 
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C. DIFFUSION WAVES.. 

' ' The bleaching wave and the thermal wave actually belong 
to a single class of diffusion waves. These diffusion waves 
appear since the rate of diffusion varies with the magnitude 
of the property that is being diffused. Specifically, diffu- 
sion proceeds more rapidly for higher values of the property 
being diffused. This process can be understood by examining 
the character of diffusion processes. 

A diffusion process is a transport process. Therefore 
two things are required? something must be transported and 
there must be a transporter. There are several examples of 
things that can be transported? energy which may be measured 
by temperature, momentum which is essentially measured by velo- 
city, and mass such as the mass of a particular species which 
is measured by concentration. There are four things that can 
perform the transporting. This first three fall into one group 
and are atoms, molecules and macroscopic particles. This 
group is simply referred to as "particles" . The last are 

massless particles, photons. 

Hence, one sees a "curious wedding of two opposite view- 
points, the continuum viewpoint and the particle viewpoint. 
Discrete particles transport — by their motion— -continuum 
properties. For example, electrons move and have a certain 
energy. The continuum viewpoint measures the average energy 
of the electrons in terms of temperature, a continuum con- 
cept. The random motion of the electrons transport energy 



214 


and the process is called thermal conduction. 

Actually, any system always has both particles and photons 
present. Thus both photon and particle transport will always 
be present? however, one or the other may dominate. 


1 . Equilibrium 

A system will always attempt to reach an equilibrium 
state. That is, the tendency is to smooth out any macroscopic 
nonuniformities that might exist. -Macroscopic nonuniformities 
are smoothed out by either macroscopic or microscopic motion 
or both. If by macroscopic motion, then a dynamical process 
ensues. If by microscopic motion, then a diffusion process 
ensues* 

Thus a diffusion process is unidirectional, it always 
proceeds toward equilibrium. This is just another statement 
of the second lav? of thermodynamics which says that on the 
average, a system of particles will always approach a more 
random disorganized state. 

a. Types of Equilibrium 

0J-0 fou?? types of equilibrium that must be estab- 
lished by diffusion or a dynamical process. The most obvious 
is spatial equilibrium of a macroscopic property. For example, 
diffusion processes attempt to smooth nonuniformities in the 
temperature . 
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The second kind is interspecies equilibrium. If there 
is more than one kind of particle in a' system, then there will 
be an exchange between the two groups of particles to establish 
equilibrium. For example a fully ionized plasma is a mixture 
of electrons and ions.’ If the electrons and ions each have a 
different temperature, then the system will seek to equalize 
the two temperatures. Another example, any system is a mix- 
ture of particles and photons. If there is an imbalance of 
energy between the photons and particles, then adjustments 
will be made to create the correct balance between the two 
systems. ‘ . • 

The third kind is intraspecies equilibrium. There is a 
certain kind of equilibrium that must be established between 
the members of a certain species. For example, suppose all 
of the members of a certain species were moving at the same 
speed. This is certainly not a equilibrium situation. Such 
a system would approach a state where all speeds are represent- 
ed. A few particles would be moving fast, and a few very slow 
and most would be somewhere in between, i.e., in equilibrium 
a "distribution" of speeds exists. An equilibrium system is 
characterized by a certain statistical distribution of its 
members. For example, there is a certain statistical distri- 
bution of energies for a system of particles which is called 
the Boltzmann distribution. This distribution function gives 
the probability of a particle having a certain energy in a 
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system with intraspacies equilibrium. -For particles in the 
classical limit it is the Boltzmann distribution that a 
particular particle species will approach within itself. 

Another examples there is a certain statistical distribution, 
function characterizing photons which is the Planck function. 

It 'gives the probability of a photon having a certain energy 
in an equilibrium system. 

The fourth kind is directional equilibrium. A system of 
particles will tend toward a state in which there is an 
equal number of particles moving in all directions with speeds 
independent of direction. This is the random component of 
velocity since macroscopic velocity always has a single direc- 
tion with a certain velocity. This equilibrium state is 
called isotropy of particle velocities. Thus four types of 
equilibrium are simultaneously being approached by the parti- 
cles in a system. Of course, the rates of approach to the 
various kinds of equilibrium are all different — but this is 
one of the things that makes diffusion processes so interesting, 
b. Approach to Equilibrium 

rf the particles of a system did not interact, then there 
would be no means by which to approach an equilibrium state; 
the particles would continue in their initial state forever. 

But particles do interact by means of collisions. There are 
four types of collisions that may occur between particles; 

1) collisions between particles in which energy or momentum 
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is shared? 2} collisions between a particle and a photon in 
which the photon is deflected (Compton scattering) ; 3) colli- 
sions between a particle and a photon in which the photon is 
extinguished (absorption)? and 4) collisions between particles 
in which a photon is emitted (emission) . 

In each type of collision , the quantity that is being 
transported is shared between the colliding particles. On the 
statistical average therefore, collisions tend to move the 
system toward equilibrium. For example? 1) collisions between 
particles tend to share energy and to promote isotropy of 
particle motion? 2) collisions between particles and photons 
promote isotropy of photon motion; 3) absorption and emission 
share particle and photon energy. 

2. Degrees of Equilibrium 

Of the four hinds of equilibrium, the degree of intra- 
species and directional equilibrium are important in discussing 
thermal and bleaching waves. 

The tendency of a group of particles is to approach a 
certain energy distribution function. A system with such an 
equilibrium distribution can be said to be in energetic equili- 
brium. The most extreme nonequilibrium situation is when all 
the particles have exactly the same energy. This state is 
called isoanergetic for particles and monochromatic for photons. 
An example of a system in energetic equilibrium is the interior 
of a perfect black body. An example of an- isoenergetic state 
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is an electron beam, and an example of a monochromatic state 
is a laser beam. 

A group of particles will approach directional equilibrium. 
Such an equilibrium is called isotropy . The extreme directional 
nonequilibrium state is when all particles are moving in exactly 
the same direction. This is called coherence . There is a 
directional state that is somewhere in between these two that 
might be called quasi— coherence . This is where all of the 
particles have a positive velocity component in say, the x 
direction, i.e., tend to be moving through the y z plane in 
the same direction. The particles are not moving in exactly 
the same direction, but in the same general direction. An • 
example of a coherent state is a laser beam^ or an electron 
beam. An example of a quasi- coherent state is the diffusion 
of particles from a very hot substance into an adjoining cold 
substance {such as the hot region behind a shock adjoining the 
cold region ahead of it). Here, the vast majority of the 
transport is crossing from the hot to the cold region. 

3. Bleaching and Thermal Waves as Diffusion Processes 

Both thermal conduction and laser irradiation are di ffu- 
sion processes. Clearly, thermal conduction is a diffusion 
process. The transporters are the particles themselves (pri- 
marily electrons in the case of a plasma) . The quantity 
being transported is energy, as measured by temperature. 
Diffusion takes place as particles migrate and collide with 
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other particles, thereby sharing their'energy. Thermal con- 
duction is usually characterized by a thermal conduction 
coefficient K and assuming that the energy flux vector 
is given by the Fourier heat conduction law, 

qr =« - K V T 

This law was derived assuming quasi-equilibrium in the energy 
both spatially and locally. A typical derivation is given 
by Vincenti and Kruger . 55 It is only valid when the spatial 
variation in the temperature is small in a mean free path, 
and when there is nearly an equilibrium distribution function. 
Thus, thermal conduction when governed by the Fourier law 
is diffusion in the equilibrium limit. 

Laser irradiation is also a diffusion process. The trans- 
'porters are photons. The quantity being transported is energy 
which may be measured by temperature, or by the number of mole- 
cules in the upper state of a two level system. Diffusion 
takes place as the photons migrate, are absorbed, and stimulate 
.the emission of other photons, thereby transporting energy into 
the system. This type of behavior which has not been identified 
as diffusion is governed by a radiative transfer equation, 

H * v - 0 


Where the absorption coefficient includes the net effect 
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of absorption and stimulated emission. Then, laser irradiation 
is diffusion in the extreme nonequilibrium limit, i.e.-, mono - 
chromatis particles moving in coherent fashion. 

Hence, laser irradiation and thermal conduction are both 
diffusion processes, which are merely at opposite ends of the 
equilibrium spectrum. Then the fact that both exhibit similar 
characteristics {bleaching waves and thermal waves) is not 
surprising. The character of equilibrium and coherent diffusion 
operating side by side is demonstrated in the next section. 

4. Rn Illustration of a Single Diffusion process 
in the Equilibrium and Coherent Limits 
Consider the example of thermal radiation in a medium as 
discussed by Zel’davich and Raiser. 58 The radiative transfer 
equation in the steady state case for a particular direction 
can be written as 


ar v 

ds 


Vv 


= K I 
Vvp 


■ £4-33) 


where X v is the spectral intensity, K v is the absorption 
coefficient (taking into account stimulated emission) . I yp 
is the spectral intensity for radiative equilibrium and s 
is the coordinate in the direction of interest. The solution 
of (4-33) is 
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X V { S ) = j“_ K v X vp exp [-f t I^ds-j ds‘ + I vo exp K v ds '] , 

(4-34) 

where I is the constant of integration., and s is the 
vo ® 

value of 3 at the edge of the system. The interpretation 
of (4-34) reveals both equilibrium and coherent diffusion. 

The first term is the radiation due to equilibrium processes. 
This term is due to incoherent radiation, which is generated 
and absorbed within the system. The second term describes the 
coherent radiation coming from outside the system. If the 
direction being considered is the same as the laser radiation, 
then. I vo is the laser intensity entering the system. In 
any other direction, i is zero. 

VO 

a quasi-coherent case would arise if the system was 
^ bounded by a much hotter black body. Then i vo would repre- ■ 
sent the spectral intensity corresponding to the temperature 
of the hot black body. x would be nonzero for all direc- 
tions that pass through the boundary into the system. 

Examination of the governing equations in both the equili- 
brium and coherent limits shows that they 'are similar in char- 
' acter even though they seemed quite different in section B. 

2n either limit, the energy equation is 



q 


{4-35} 
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In the coherent limit, the intensity is given by 


I = I Q exp c-r K ds’J (4-36) 

s Q 

The energy flux is the radiation intensity; 3 = l e g where 

e is a unit vector in the direction of the laser beam, 
s 

In the equilibrium limit, the energy equation is the 
same as (4-35) but the energy flux vector is different. The 
spectral energy flux vector S v is given by 






vu 


vp 


(4-37) 


where is the absorption length at frequency v , and 

u vp tJle e | 3 u il^ 1:>r i um spectral radiant energy density (the 
Planck function). This equation, called the diffusion approx- 
imation comes from assuming quasi equilibrium locally and 
spatially. 5 ? But the total energy flux is needed and this 
requires an integral of (4-37) over all frequencies v . The 
integration yields the energy flux for the equilibrium limit. 


S 




VT 


(4-38) 


/ is the Rosseland mean free path, an averaged value of L v 

(with proper weighting). 58 U p is the radiant energy density 

and is given by U = 4aT 4 /c where a is the Stefan Boltzmann 
P 
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constant. So then, in the equilibrium limit, the radiation 
and absorption mechanism actually behaves like thermal energy, 
diffusing" with a conduction coefficient, 

K = ^ 1 0 T 3 (4-39) 

This radiation process has been identified as a diffusion 
phenomenon and shown to result in thermal waves by Zel’dovich 

„ .. 54 

and Raizer. 

i 

The similarity of the energy flux in the coherent (4-34) 
ana the equilibrium (4-38) limits is easily seen* In the 
equilibrium limit consider a hot body bordering a cold body 
at x = 0 . The radiation from the hot body diffuses into 
the cold body. As the temperature of the cold body increases, 
the coefficient of radiation thermal conduction (proportional 
to T 3 as seen in (4—39) ) increases greatly- So a thermal 

wave phenomenon will occur similar to that in section B. 

In the coherent limit, consider a laser beam which enters 
the cold body at x = 0. If the absorption coefficient K v 
decreases with increasing temperature (as with inverse brem- 
strahlung or any other bleachable type absorption) , then the 
bleaching wave phenomenon will arise. 

Thus, these two limits of the same physical phenomenon 
produce similar results. Granted, the details of the processes 
may vary somewhat but the general character is the same. The 
next section discusses a number of physical phenomena that 
are governed by nonlinear diffusion. Among them are examples 
in both the equilibrium and coherent limits. 
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D. EXAMPLES OF DIFFUSION WAVE PHENOMENA 
There are many physical examples of nonlinear diffusion 
phenomena, and almost without exception the nonlinearity is 
'of the bleachahle variety so that diffusion waves appear. 
Several physical examples are discussed below. A chart sura 
marizing these phenomena is given in figure 61. 

1. Particle Diffusion 
a. Diffusion of Mass or Species 
The diffusion of a plasma through a magnetic field is 
discussed by Kaufman. 59 The property being diffused is mass 
as measured by the number density of the particles. The 
diffusion coefficient is proportional to the number density, 
so this is bleachable diffusion . Kaufman presents a similar- 
ity solution that clearly demonstrates the diffusion wave. 
This diffusion is the "equilibrium" variety since it is char- 
acterized by a diffusion coefficient. 

The diffusion of lead chloride in solid silver chloride 
is studied by Wagner. 60 The quantity being diffused is the 
molecules of a certain species and is measured as concentra- 
tion. The diffusion coefficient is proportional to the 
concentration, so then it is governed by the same equation as 
the diffusion of a plasma in a magnetic field. 

The diffusion of H 2 0 in swelling membranes was studied 
by Fatt and Goldsti'ck. 61 In. this instance the diffusion 
coefficient decreases with increasing concentration so that 
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Figure 61. Examples of diffusion phenomena 
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the diffusion is not the bleachable variety. Then no diffu- 
sion wave appears. This was the only instance found of a 
nonbleachable diffusion mechanism although there are pro- 
bably ether examples. This is diffusion in the equilibrium 
limit since a diffusion coefficient is used. 

b. Diffusion of Energy 

The diffusion of thermal energy by nonlinear thermal 
conduction has heen demonstrated in the numerical studies of 

CO 

Shearer, and in the analytic studies of Babuel-Peyrissae, 

CO 

et. al. The phenomenon arose with the rapid laser heating 
of an overdense plasma. The laser radiation only penetrates 
to the critical density, hence only heats the underdense 
region directly. But the hot underdense region acts like an 
energy source and rapidly conducts energy to the overdense 
“region. The thermal conduction coefficient is proportional 
to the five halves power of temperature, so is bleachable. 
Thus, a diffusion wave moves into the overdens.e material. 
Interestingly enough, in the calculation of Babuel-Peyrissae, 
et. al., the thermal wave moved very fast at first then slowed 
down, and was succeeded by a shock wave propagating into the 
material. This illustrates that diffusion waves can move 
faster than dynamical processes in, some cases. 

The preheating by thermal conduction in front of a strong 

54 

shock wave is discussed by Zel'dovich and Raizer. This is 
a steady state diffusion phenomena so the temperature profiles 
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produces are stationary relative to the shock wave. Never- 
theless , the familiar humped shape in'-the temperature profile 
arises due to the bleachable diffusion. If the thermal con- 
duction were constant, then only the exponentially decaying 
profile would he observed. 

'Heating by an electron beam is suggested to be an example 
of coherent diffusion. This is conceptually the same as the 
heating by-a coherent monochromatic beam of photons, i.e., 
a laser. Some details may differ since electron beams generate 
streaming effects leading to the two stream instability, some- 
thing that may not have an analogy in photon beams. 

2. Photon Diffusion: ' Diffusion of Energy 

The laser heating of a stationary plasma has been studied 
in chapter II. This is an example of coherent diffusion. The 
diffusion is bleachable so the familiar bleaching wave appears. 
Laser irradiation of a two level substance has also been 
discussed in section A 1. It is obviously a bleachable 
diffusion phenomenon so the bleaching wave appears. Any other 
bleachable absorption mechanism, such as bound free absorption, 
will also result in a bleaching wave. 

The ^breakdown wave, an example of a backwards moving 
bleaching wave (or an unbleaching wave) was also discussed 
in section A 2. This is also coherent diffusion. The reverse 
nature of the wave is due to the converging geometry of a 
focused laser beam-. 
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The initial expanding fireball in an atomic explosion 
has been explained by Marshak 64 as the diffusion of equili- 
brium radiation. In this example, the radiation diffuses 
rapidly at first, then slows down and is succeeded by a shock 
wave— in the same manner as in the study of Babuel-Peyrissae, 
et. al. However in this instance it is equilibrium radiation 
that is diffusing. So this is an example of equilibrium 
diffusion by photons. The equations governing this problem 
were mentioned in section C 4: (4-33) and (4-36). 

These are a few examples of nonlinear diffusion ana 
doubtless there are others. But it is clear that many 
physical phenomena— especially those characterized by high 
temperatures— are characterized by the diffusion wave. 
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